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Abstract 

We consider generalized Wigner ensembles and general /3-ensembles with analytic potentials 
for any P > 1. The recent universality results in particular assert that the local averages of 
consecutive eigenvalue gaps in the bulk of the spectrum are universal in the sense that they 
coincide with those of the corresponding Gaussian /3-ensembles. In this article, we show that 
local averaging is not necessary for this result, i.e. we prove that the single gap distributions 
in the bulk are universal. 
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1 Introduction 



The fundamental vision that random matrices can be used as basic models for large quantum 
systems was due to E. Wigner [53]. He conjectured that the eigenvalue gap distributions of large 
random matrices were universal ( "Wigner surmise" ) in the sense that large quantum systems and 
random matrices share the same gap distribution functions. The subsequent works of Dyson, 
Gaudin and Mehta clarified many related issues regarding this assertion and a thorough under- 
standing of the Gaussian ensembles has thus emerged (see the classical book of Mehta [40] for a 
summary). There are two main categories of random matrices: the invariant and the non- invariant 
ensembles. The universality conjecture, which is also known as the Wigner-Dyson-Gaudin-Mehta 
(WDGM) conjecture, asserts that for both ensembles the eigenvalue gap distributions are universal 
up to symmetry classes. For invariant ensembles, the joint distribution function of the eigenvalues 
can be expressed explicitly in terms of one dimensional particle systems with logarithmic interac- 
tions (i.e., log-gases) at an inverse temperature /3. The values /3 = 1,2,4, correspond to the classical 
orthogonal, unitary and symplectic ensembles, respectively. Under various conditions on the ex- 
ternal potential, the universality for the classical values (3 = 1,2,4 was proved, via analysis on the 
corresponding orthogonal polynomials, by Fokas-Its-Kitaev [31], Deift et. al. [11, 14, 15], Bleher- 
Its [6], Pastur-Shcherbina [43, 44] and in many consecutive works, see e.g. [12, 13, 38, 45, 52]. For 
nonclassical values of /3 there is no matrix ensemble behind the model, except for the Gaussian 
cases [17] via tridiagonal matrices. One may still be interested in the local correlation functions of 
the log-gas as an interacting particle system. The orthogonal polynomial method is not applicable 
for nonclassical values of /3 even for the Gaussian case. Nevertheless, the local statistics can still be 
described very precisely in the Gaussian case with a different method [50, 51]. The universality for 
general /3-ensembles was established only very recently [7, 8] by a new method based on dynamical 
methods using Dirichlet form estimates from [21, 22]. This method is important for this article and 
we will discuss it in more details later on. All previous results achieved by this method, however, 
required in their statement to consider a local average of consecutive gaps. In the current paper 
we will prove universality of each single gap in the bulk. 

Turning to the non-invariant ensembles, the most important class is the N x N Wigner matrices 
characterized by the independence of their entries. In general, there is no longer an explicit 
expression for the joint distribution function for the eigenvalues. However, there is a special 
class of ensembles, the Gaussian divisible ensembles, that interpolate between the general Wigner 
ensembles and the Gaussian ones. For these ensembles, at least in the special Hermitian case, 
there is still an explicit formula for the joint distribution of the eigenvalues based upon the Harish- 
Chandra-Itzykson-Zuber integral. This formula was first put into a mathematically useful form 
by Johansson [36] (see also the later work of Ben Arous-Pechc [5]) to prove the universality of 
Gaussian divisible ensembles with a Gaussian component of size order one. In [19], the size of 
the Gaussian component needed for proving the universality was greatly reduced to N~^^'^~^^. 
More importantly, the idea of approximating Wigner ensembles by Gaussian divisible ones was 
first introduced and, after a perturbation argument, this resulted in the first proof of universality 
for Hermitian ensembles with general smooth distributions for matrix elements. The smoothness 
condition was later on removed in [46, 20]. 

In his seminal paper [18], Dyson observed that the eigenvalue distribution of Gaussian divisible 
ensembles is the same as the solution of a special system of stochastic differential equations, 
commonly known now as the Dyson Brownian motion, at a fixed time t. For short times, t is 
comparable with the variance of the Gaussian component. He also conjectured that the time to 
"local equilibrium" of the Dyson Brownian motion is of order 1/iV, which is then equivalent to the 
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universality of Gaussian divisible ensembles with a Gaussian component of order slightly larger 
than N-^/"^. Thus the work [19] can be viewed as proving Dyson's conjecture for the Hermitian 
case. This method, however, completely tied with an explicit formula that is so far restricted to 
the Hermitian case. 

A completely analytic approach to estimate the time to local equilibrium of the Dyson's Brow- 
nian motion was initiated in [21] and further developed in [22, 26, 25], see [30] for a detailed 
account. In these papers, Dyson's conjecture in full generality was proved [26] and universality 
was established for generalized Wigner ensembles for all symmetric classes. The idea of a dynami- 
cal approach in proving universality turns out to be a very powerful one. Dyson Brownian motion 
can be viewed as the natural gradient flow for Gaussian /3 log-gases (we will often use the termi- 
nology /3 log-gases for the /3-ensembles to emphasize the logarithmic interaction). The gradient 
flow can be defined with respect to all /? log-gases, not just the Gaussian ones. Furthermore, one 
can consider gradient flows of local log-gases with fixed "good boundary conditions" . Here "local" 
refers to Gibbs measures on 7^,0 < a < 1, consecutive points of a log-gas with the locations of 
all other points fixed. By "good boundary conditions" we mean that these external points are 
rigid^ i.e. their locations arc close to their classical locations given by the limiting density of the 
original log-gas. Using this idea, we have proved the universality of general /3-cnscmbles in [7, 8] 
for analytic potential. 

The main conclusion of these works is that the local gap distributions of cither the generalized 
Wigner ensembles (in all symmetry classes) or the general /3-ensembles are universal in the bulk of 
the spectrum (see [29] for a recent review). The dynamical approach based on Dyson's Brownian 
motion and related flows also provides a conceptual understanding for the origin of the universality. 
For technical reasons, however, these proofs apply to averages of consecutive gaps, i.e. cumulative 
statistics of A^^ consecutive gaps were proven to be universal. Averaging the statistics of the 
consecutive gaps is equivalent to averaging the energy parameter in the correlation functions. 
Thus, mathematically, the results were also formulated in terms of universality of the correlation 
functions with averaging in an energy window of size A^"^'*''^. 

The main goal of this paper is to remove the local averaging in the statistics of consecutive 
gaps in our general approach using Dyson Brownian motion for both invariant and non-invariant 
ensembles. Wc will show that the distribution of each single gap in the bulk is universal, which 
we will refer to as the single gap universality or simply the gap universality whenever there is no 
confusion. The single gap universality was proved for a special class of Hermitian Wigner matrices 
with the property that the first four moments of the matrix elements match those of the standard 
Gaussian random variable [47] and no other results have been known before. In particular, the 
single gap universality has not been proved even for the Gaussian orthogonal ensemble (GOE). 

The gap distributions are closely related to the correlation functions which were often used to 
state the universality of random matrices. These two concepts are equivalent in a certain average 
sense. However, there is no rigorous relation between correlation functions at a fixed energy and 
single gap distributions. Thus our results on single gap statistics do not automatically imply the 
universality of the correlation functions at a fixed energy which was rigorously proved only for 
Hermitian Wigner matrices [19, 46, 20, 30]. 

The removal of a local average in the universality results proved in [21, 27, 28] is a technical 
improvement in itself and its physical meaning is not especially profound. Our motivation for 
taking seriously this endeavor is due to that the single gap distribution may be closely related to 
the distribution of a single eigenvalue in the bulk of the spectrum [33] or at the edge [48, 49]. Since 
our approach does not rely on any explicit formula involving Gaussian matrices, some extension 
of this method may provide a way to understand the distribution of an individual eigenvalue of 
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Wigner matrices. 

The main new idea in this paper is an analysis of the Dyson Brownian motion via parabohc 
regularity using the De Giorgi-Nash-Moser idea. Since the Hamiltonians of the local log-gases 
are convex, the correlation functions can be re-expressed in terms of a time average of certain 
random walks in random environments. The connection between correlation functions of general 
log-concave measures and random walks in random environments was already pointed out in the 
work of Hclffer and Sjostrand [35] and Naddaf and Spencer [41]. This connection was used as an 
effective way to estimate correlation functions for several models in statistical physics, see, e.g. 
[2, 16, 32, 34]. 

In this paper we observe that the single gap universality is a consequence of the Holder regularity 
of the solutions to these random walk problems. Due to the logarithmic interaction, the random 
walks are long ranged and their rates may be singular. Furthermore, the random environments 
themselves depend on the gap distributions, which were exactly the problems we want to analyze! If 
we view these random walks as (discrete) parabolic equations with random coefficients, we find that 
they are of divergence form and are in the form of the equations studied in the fundamental paper 
by Caffarelli, Chan and Vasseur [10]. The main difficulty to apply [10] to gain partial regularity is 
that the diffusion coefficients in our settings are random and can be much more singular than were 
allowed in [10]. For controlling the singularities of these coefficients, we prove an optimal level 
repulsion estimate for the local log-gases. With these estimates, we are able to extend the method 
of [10] to prove Holder regularity for the solution to these random walks problems. This shows 
that the single gap distributions are universal for local log-gases with good boundary conditions, 
which is the key result of this paper. 

For /3-ensembles, it is known that the rigidity of the eigenvalues ensures that boundary con- 
ditions are good with high probability. Thus we can apply the local universality of single gap 
distribution to get the single gap universality of the /3-ensembles. We remark, however, that the 
current result holds only for /3 > 1 in contrast to /3 > in [7, 8], since the current proof heavily re- 
lies on the dynamics of the gradient flow of local log-gases. For non-invariant ensembles, a slightly 
longer argument using the local relaxation flow is needed to connect the local universality result 
with that for the original Wigner ensemble. This will be explained in Section 5. 

In summary, we have recast the question of the single gap universality for random matrices, 
envisioned by Wigner in the sixties, into a problem concerning the partial regularity of a parabolic 
equation in divergence form studied by De Giorgi-Nash-Moser. Thanks to the insight of Dyson 
and the important progress by Caffarelli-Chan- Vasseur [10], we are able to establish the WDGM 
universality conjecture for each individual gap via De Giorgi-Nash-Moser's idea. We now introduce 
our models rigorously and state the main results. 

2 Main results 

We will have two related results, one concerns the generalized Wigner ensembles, the other one the 
general /3-ensembles. We first define the generalized Wigner ensembles. Let H = {hij)f^j^i be an 
N X N hermitian or symmetric matrix where the matrix elements hij = hji, i < j, are independent 
random variables given by a probability measure Uij with mean zero and variance afj > 0; 

E/iy=0, 4 E|/i,,f . (2.1) 
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The distribution Vij and its variance af^ may depend on N , but we omit this fact in the notation. 
We also assume that the normahzed matrix elements have a uniform subexponential decay, 

> x<T.,j) < 6'icxp(-x*^), X > 0, (2.2) 

with some fixed constants 9i,92 > 0, uniformly in N,i, j. In fact, with minor modifications of the 
proof, an algebraic decay 

P(|/iy| > xa^,) < Cmx~^' 
with a large enough M is also sufficient. 

Definition 2.1 ([24]) The matrix ensemble H defined above is called generalized Wigner matrix 
if the following assumptions hold on the variances of the matrix elements (2.1) 



(A) For any j fixed 



N 



E4=l- (2-3) 



i=i 



(B) There exist two positive constants, Cm/ and Csup, independent of N such that 

N - - N ^ ' 

Let P and E denote the probability and the expectation with respect to this ensemble. 

We will denote by Ai < A2 < . . . < Aat the eigenvalues of H. In the special case when afj = 1/N 
and hij is Gaussian, the joint probability distribution of the eigenvalues is given 

-NI3U(\) ^ \2 1 

A. = M[f ^(dA) = dA, n{\) = ^ ^ - - ^ log |A, - A,;|. (2.5) 

The value of /3 depends on the symmetry class of the matrix; (3 = 1 for GOE, (3 = 2 for GUE and 
13 = A for GSE. Here Zp is the normalization factor so that /i is a probability measure. 

It is well known that the density or the one point correlation function of /i converges, as N 00, 
to the Wigner semicircle law 

q{x) := -^V(4-2;2)+. (2.6) 

ZTT 

We use the notation 7^ for the j-th quantile of this density, i.e. 7j is defined by 

QG{x)dx. (2.7) 

2 



]_ 

N 



Our main result on the generalized Wigner matrices asserts that the local gap statistics in the 
bulk of the spectrum are universal for any general Wigner matrix, in particular they coincide with 
those of the Gaussian case. 

For any integers A < B we introduce the notation [A, B] := {A, A + 1, . . . , _B}. 

Theorem 2.2 (Gap universality for Wigner matrices) Let H be a generalized Wigner en- 
semble with subexponentially decaying matrix elements, (2.2). Fix a positive number a > 0, an 



6 



integer n G N and a smooth, compactly supported function O : M" — > M. There exists an e > and 
C > 0, depending only on a and on O such that 

[E - E^]0(7V(xj - Xj+i),N{xj - Xj+2), . . .,N{xj - Xj+n)) \ < CN-', (2.8) 

for any j € faN, (1 — a)-/V]l and for any sufficiently large N > Nq, where Nq depends on all 
parameters of the model, as well as on n and a. 

More generally, for any fc, m G faN, (1 — a)-^] we have 

EO{iNQk)ixk - Xk+i), {Ngk){xk - Xk+2), {Ngk){xk - Xk+„)) (2.9) 



),..., {Ng,n){Xm - Xm+n)) 

where the local density g^ is defined by :~ gijk)- 



< CN- 



It is well-known that the gap distribution of Gaussian random matriees for all symmetry classes 
can be explicitly expressed via a Fredholm determinant provided that a certain local average is 
taken, see [11, 12, 13]. The result for a single gap, i.e. without local averaging, was only achieved 
recently in the special case of the Gaussian unitary ensemble (GUE) by Tao [47] (which then easily 
implies the same results for Hermitian Wigner matrices satisfying the four moment matching 
condition). It is not clear if a similar argument can be applied to the GOE case. 

We now define the /3-ensembles with a general external potential. Let /3 > be a fixed param- 
eter. Let V{x) be a real analytic potential on R that grows faster than (2 + e) log |x| at infinity 
and satisfies 

infy">-oo. (2.10) 

R 

Consider the measure 



= /i(^^(dA) = ^ dA, n{X) = -Y,V{X,) - -Y,log\\, ^ K\. (2.11) 

i=l i<j 



Since fj, is symmetric in all its variables, we will mostly view it as a measure restricted to the 
simplex 

S(^) {A : Ai < A2 < . . . < Ajv} C M^. (2.12) 

Note that the Gaussian measure (2.5) is a special case of (2.11) with V{X) ~ A^/2. In this case we 
use the notation nc for /z. 
Let 

1 ^ 

gW(A) :=E''-5]<5(A-A,) 

denote the density, or the one-point function, of//. It is well known [1, 9] that g^^'' converges weakly 
to the equilibrium density g = gv as iV — >■ 00. The equilibrium density can be characterized as 
the unique minimizer (in the set of probability measures on R endowed with the weak topology) 
of the functional 

/(zy) = / V{t)diy{t)- II \og\t- s\di^{t)di'{s). (2.13) 
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In the case, V{x) = the mimmizer is the Wigner semicircle law g ~ go, defined in (2.6), 

where the subscript G refers to the Gaussian case. In the general case we assume that p = is 
supported on a single compact interval, [A,B] and g G C'^{A,B). Moreover, we assume that V 
is regular in the sense that g is strictly positive on [A, B) and vanishes as a square root at the 
endpoints, see (1.4) of [8]. It is known that these condition are satisfied if, for example, V is strictly 
convex. 

For any j < N define the classical location of the j-th particle 'jjy by 

= gvix)dx, (2.14) 

and for the Gaussian case we have [A, B] = [—2, 2] and we use the notation jj,G = Ij for the 
corresponding classical location, defined in (2.7). We set 

qJ Qvilj.y), and gf Qoilj-G) (2.15) 

to be the limiting density at the classical location of the j-th particle. Our main theorem on the 
^-ensembles is the following. 

Theorem 2.3 (Gap universality for /3-ensembles) Let /3 > 1 and V be a real analytic poten- 
tial with (2.10) such that gv is supported on a single compact interval, [A,B], gv € C^(A, _B), 
and that V is regular. Fix a positive number a > 0, an integer n g N and a smooth, compactly 
supported function O : M" — >■ M. Let /i ~ fiy ~ l^\fv given by (2.11) and let fiQ denote the 
same measure for the Gaussian case. Then there exist an e > 0, depending only on a, /3 and the 
potential V , and a constant C depending on O such that 

W^-0(^iNgl)ixk - Xk+i), {NgDixk - Xk+2), ■ • • , (NgD^Xk - Xk+n)) (2.16) 

- E'"='0(^{Ng^){x^ - x„^+l), {Ng^){x^ - x^+s), . . . , iNgg){x,n ~ x„^+n)) < 

for any k,m ^ \olN, (1 — oi)N\ and for any sufficiently large N > Nq, where Nq depends on V, 
(3, as well as on n and a. In particular, the distribution of the rescaled gaps w.r.t. fiy does not 
depend on the index k in the bulk. 

Theorem 2.3, in particular, asserts that the single gap distribution in the bulk is independent 
of the index k. The special GUE case of this assertion is the content of [47] where the proof uses 
some special structures of GUE. 

The proofs of both Theorems 2.2 and 2.3 rely on the uniqueness of the gap distribution for a 
localized version of the equilibrium measure (2.5) with a certain class of boundary conditions. This 
main technical result will be formulated in Theorem 3.1 in the next section after we introduce the 
necessary notations. An orientation of the content of the paper will be given at the end of Section 
3.1. 

We remark that Theorem 2.3 is stated only for /3 > 1; on the contrary, the universality with 
local averaging in [7, 8] was proved for /3 > 0. The main reason is that the current proof relies 
heavily on the dynamics of the gradient flow of local log-gases. Hence the well-posedness of the 
dynamics is crucial which is available only for /3 > 1. On the other hand, in [7, 8] we use only 
certain Dirichlet form inequalities (see, e.g. Lemma 5.9 in [7]), which we could prove with an 
effective regularization scheme for all /? > 0. For /3 < 1 it is not clear if such a regularization can 
also be applied to the new inequalities we will prove here. 
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3 Local equilibrium measures 

3.1 Basic properties of local equilibrium measures 

Fix two small positive numbers, a,6 > 0. Choose two positive integer parameters L, K such that 

L e laN, (1 - a)iV|, < K < TV^/". (3.1) 

We consider the parameters L and K fixed and often we will not indicate them in the notation. 
All results will hold for any sufficiently small a, 5 and for any sufficiently large N > Nq, where 
the threshold A^o depends on a, S and maybe on other parameters of the model. Throughout the 
paper we will use C and c to denote positive constants which, among others, may depend on a, S 
and on the constants in (2.2) and (2.4), but we will not emphasize this dependence. Typically C 
denotes a large generic constant, while c denotes a small one whose values may change from line to 
line. These constants are independent of K and N, which are the limiting large parameters of the 
problem, but they may depend on each other. In most cases this interdependence is harmless since 
it only requires that a fresh constant C be sufficiently large or c be sufficiently small, depending on 
the size of the previously established generic constants. In some cases, however, the constants are 
related in a more subtle manner. In this case we will use Co, Ci, . . . and cq, ci, . . . etc. to denote 
specific constants in order to be able to refer to them along the proof. 
For convenience, we set 

K. := 2/^ + 1. 

Denote / = Il k '■= l^ — K, L + K^ the set of K. consecutive indices in the bulk. We will distinguish 
the inside and outside particles by renaming them as 

(Ai,A2,...,AAr) := {yi, . . .yL^K~i,XL^K, ■ ■ ■ ,XL+K,yL+K+i, ■ ■ - Vn) e S"^^'. (3.2) 

Note that the particles keep their original indices. The notation S'^' refers to the simplex (2.12). 
In short we will write 

{xl~k,---xl+k), and y = {yi, ■ ■ ■yL-K-i,yL+K+i, ■ ■ -Vn)- 

These points are always listed in increasing order, i.e. x G S''^' and y G We will refer to 

the y's as the external points and to the x's as internal points. 

We will fix the external points (often called boundary conditions) and study the conditional 
measures on the internal points. We first define the local equilibrium measure (or local measure in 
short) on x with boundary condition y by 

(3.3) 

where fj, = /^(y,x) is the (global) equilibrium measure (2.11) (we do not distinguish between the 
measure /i and its density function ^(y,x) in the notation). Note that for any fixed y S ~_iN-K.)^ 
all Xj lie in the open configuration interval, denoted by 

J = Jy := iyL-K-i,yL+K+i)- 
Define ^ 

y ■■= -^{yL-K-i + VL+K+i) 
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/^y(dx) %(x)dx, ^y(x) := M(y,x) 



j M(y,x)dx 



to be the midpoint of the configuration interval. We also introduce 

aj:=y+^^l^l, 3&Il^k, (3.4) 

to denote the /C equidistant points within the interval J. 

For any fixed L, K, y, the equilibrium measure can also be written as a Gibbs measure, 

My = = (3-5) 

with Hamiltonian 



Vy{x):=V{x)~^J2^og\x-y,\. (3.6) 

Here V^(a;) can be viewed as the external potential of a /3- log-gas of the points {xi : i e /} in the 
configuration interval J. 

Our main technical result. Theorem 3.1 below, asserts that, for K,L chosen according to (3.1), 
the local gap statistics is essentially independent of V and y as long as the boundary conditions y 
are regular. This property is expressed by defining the following set of "good" boundary conditions 
with some given positive parameters v, a: 

■Rlm = nL,K{^, a) :={y : \yk - 7fc| < ^"'+", k e {aN, (1 - a)Nl \ Il,k} (3.7) 

n {y : \yk - 7^1 < N-^/^''+'', k e [7v3/5+.^ ^ „ iv3/5+-j} 
n {y : \yk ' 7k\ < 1, k e U^NjX Il,k}- 

In Section 4 we will see that this definition is taylored to the previously proven rigidity bounds for 
the /3-ensemble, see (4.4). The rigidity bounds for the generalized Wigner matrices are stronger, 
see (5.1), so this definition will suit the needs of both proofs. 

Theorem 3.1 (Gap universality for local measures) Fix L,L and JC = 2K + 1 satisfying 
(3.1) with an exponent S > 0. Consider two boundary conditions y,y such that the configuration 
intervals coincide, 

J = iyL-K-i,yL+K+i) = iyi-K-i^vi+K+i)- (3-8) 

We consider the measures fi ~ /iy^^.y and Jl ~ fi~ ^ y defined as in (3.5), with possibly two different 
external potentials V and V . Let ^ > 6e a small constant. Assume that \J\ satisfies 



Suppose that y,y £ 7^l.x(C^'5/2, a/2) and that 



max 



E^'yxJ-aJ + max E^^x^ - <CN-^K^ (3.10) 
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holds. Let the integer number p satisfy \p\ < K — K^^^ for some small ^* > 0. Then there 
exists ^0 > 0, depending on S, such that if S,,^* < then for any n fixed and any bounded smooth 
observable O : M" — >■ M with compact support we have 



E''yO{N{xL+p - XL+p+l), ■ . -N^XL+p - XL+p+n)) 



(3.11) 



yO{N{xj^ - xj^ ), . . . N{x~ - XI )) 



< CK 



for some e > depending on 5, a and for some C depending on O. This holds for any N > Nq 
sufficiently large, where Nq depends on the parameters and C in (3.10). 

In the following two theorems we establish rigidity and level repulsion estimates for the local 
log-gas fiy with good boundary conditions y. While both rigidity and level repulsion are basic 
questions for log gases, our main motivation to prove these theorems is to use them in the proof of 
Theorem 3.1. The current form of the level repulsion estimate is new which has not been proved in 
literature (a much weaker form was proved in (4.11) of [7]). The rigidity estimate was proved for 
the global equilibrium measure /i in [7]. From this estimate, one can conclude that fj,y has a good 
rigidity bound for a set of boundary conditions with high probability w.r.t. the global measure 
fi. However, we will need a rigidity estimate for fiy for a set of y's with high probability with 
respect to some different measure, which may be asymptotically singular to /j, for large N. For 
example, in the proof for the gap universality of Wigner matrices such a measure is given by the 
time evolved measure ftfi, see Section 5. The following result asserts that a rigidity estimate holds 
for /iy provided that y itself satisfies a rigidity bound and an extra condition, (3.12), holds. This 
provides explicit criteria to describe the set of "good" y's whose measure w.r.t. /t/i can then be 
estimated with different methods. 

Theorem 3.2 (Rigidity estimate for local measures) Let L and K satisfy (3.1) with 6 the 
exponent appearing in (3.1). Let ^, a be any fixed positive constants. For y G TZLji{^S/2,a) 
consider the local equilibrium measure /iy defined in (3.5) and assume that 



(3.12) 



is satisfied. Then there are positive constants C, c, depending on ^, such that for any k ^ I and 
u>0, 

P^^(|a;fe - ak\ > uK^N'^^ < Ce~™'. (3.13) 
Now wc state the level repulsion estimates which will be proven in Section 6.2. 

Theorem 3.3 (Level repulsion estimate for local measures) Let L and K satisfy (3.1) and 
let ^, a be any fixed positive constants. Then for y G TZl.k ~ Ti-L.K{S,'^6 /2, a) we have the following 
estimates: 

i) [Weak form of level repulsion] For any s > we have 



and 



[x,+i -x^< s/N]<C {Nsf^^ , 
^^-[x,+2 ~x,< s/N] < C{Nsf''+^ 



ielL-K-l,L + K} 
ielL-K-l,L + K-l\ 



(3.14) 
(3.15) 
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(Here we used the convention that xl-k-i '■= J/i-X-i, •'J^L+A'+i := Vl+k+i-) 

a) [Strong form of level repulsion] Suppose that there exist positive constants C,c such that the 
following rigidity estimate holds for any k Cz I : 

P^^(|a;fc-afc| > Cif«'iV-i) <Ccxp(-i^"). (3.16) 

Then there exists small a constant 9, depending on C, c in (3.16), such that for any s > cxp(— A"^). 
we have 

P''^[x,+i-x,<s/N]<C{K^s\ogNf^\ ielL-K -1,L + Kl (3.17) 

and 

¥''y[x^+2-x, <s/N]<C {K^s log Nf^^\ i E IL - K - 1, L + K - Ij. (3.18) 



We remark that the estimates (3.18) and (3.15) on the second gap are not needed for the main 
proof, we hsted them only for possible further reference. The exponents are not optimal; one would 
expect them to be 3/3 + 3. With some extra work, it should not be difficult to get the optimal 
exponents. Moreover, our results can be extended to Xi+k — Xi for any k fin We also mention that 
the assumption (3.16) required in part ii) is weaker than what we prove in (3.13). In fact, the 
weaker form (3.16) of the rigidity would be enough throughout the paper except at one place, at 
the end of the proof of Lemma 6.8. 

Theorem 3.1 is our key result. In Sections 4 and 5 we will show how to use Theorem 3.1 to 
prove the main Theorems 2.2 and 2.3. Although the basic structure of the proof of Theorem 2.3 is 
similar to the one given in [7] where a locally averaged version of this theorem was proved under 
a locally averaged version of Theorem 3.1, here we have to verify the assumption (3.10) which will 
be done in Lemma 4.1. The proof of Theorem 2.2, on the other hand, is very different from the 
recent proof of universality in [27, 28]. This will be explained in Section 5. 

The proofs of the auxiliary Theorems 3.2 and 3.3 will be given in Section 6. The proof of 
Theorem 3.1 will start from Section 6.3 and will continue until the end of the paper. At the 
beginning of Section 6.3 we will explain the main ideas of the proof. For readers interested in the 
proof of Theorem 3.1, Sections 4 and 5 can be skipped. 



3.2 Extensions and further results 

We formulated Theorems 3.1, 3.2 and 3.3 with assumptions requiring that the boundary conditions 
y are "good" . In fact, all these results hold in a more general setting. 

Definition 3.4 An external potential U of a j3-log-gas of K points in a configuration interval 
J = (a, h) is called -regular, if the following hounds hold: 

C/'(.T) = f?(y)log^ + of-^), xeJ, (3.20) 



d-{x) \Nd{x) 

U"{x) > inf + x£j, (3.21) 

a(x) 
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with some positive c > and for some small ^ > 0, where 

d{x) := niin{|a; — a|, |a; — 6|} 

is the distance to the boundary of J and 

d_{x) := d{x) + g{y)N-^Ki, d+{x) := niax{|a- - a\, \x - b\} + g{y)N-'^K^. 

The following lemma, proven in Appendix A, asserts that "good" boundary conditions y give 
rise to regular external potential Vy. 

Lemma 3.5 Let L and K satisfy (3.1) and S is the exponent appearing in (3.1). Then for any 
y G 7?.i_iv'(C^/2, Q;/2) the external potential Vy (3.6) on the configuration interval Jy is -regular; 

F;(x) = ,(S)l„g^±g>+0(-^). XCJ„ (3^23) 

V^{x) > iniV" + x e Jy. (3.24) 

Inspecting the proofs of Theorems 3.1, 3.2, 6.3 and Theorem 3.3, they use only the property 
that Vy on Jy is regular. Thus the same proofs are valid for any local /3-log-gas with an external 
potential U, where the conditions y G TZL.K{^S/2,a) or y £ TIl.k{^'^S /2, a) in these statements 
are replaced by requiring that U is A'*-regular or A'* -regular, respectively. 



4 Gap universality for /3-ensembles: proof of Theorem 2.3 
4.1 Rigidity bounds and its consequences 

The aim of this section is to use Theorem 3.1 to prove Theorem 2.3. In order to verify the 
assumptions of Theorem 3.1, we first recall the rigidity estimate w.r.t. fj, defined in (2.11). Recall 
that 7fc = jk,v denotes the classical location of the fc-th point (2.14). For the case of convex 
potential, in Theorem 3.1 of [7] it was proved that for any fixed a > and i' > 0, there are 
constants Co, ci, C2 > such that for any A^ > 1 and k G laN, (1 — a)A^l, 

P^{\Xk-lk\> N-^+'') <Cocxp{-ciN--'). (4.1) 

The same estimate holds also for the non-convex case, see Theorem 1.1 of [8]. 

Near the spectral edges, a somewhat weaker control was proven, see Lemma 3.6 of [7]: for any 
J/ > there are Co , ci , C2 > such that 

P'^ {\Xk-lk\ > A^-4/i5+-) < Coexp(-ciAr^-) (4.2) 

for any N^/^+'' <k<N- N^/^+'', if N > Noii^) is sufficiently large. We can choose Co, ci, C2 to 
be the same in (4.1) and (4.2). 

We claim that the estimate (4.2) holds also for the non-convex case. Inspecting the proof of 
Lemma 3.6 of [7], the convexity of V was used only in the logarithmic Sobolev inequality (LSI) 
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three times. The first apphcation of LSI was in the anaiysis of the loop equation. As it was already 
observed in [8], Lemma 2.1 of [8] can replace LSI for this purpose, in particular Lemma 2.2 of [8] 
holds without convexity. The second application of LSI was to obtain a concentration estimate 
on \j — K^Xj on scale iV~^/^ uniformly in j (see (3.12) of [7]). In Lemma 3.6 of [8], however, 
the equivalence of the original measure ^ and its convexified version i/ (see Definition 3.4 of [8]) 
was established. In particular, the same concentration estimate was established for the non-convex 
measure as well (see (3.22) of [8]) and the analogue of Corollary 3.5 of [7] was also obtained (see 
(3.24) of [8]). Finally, LSI was used also to control the variance term (3.32) of [8]. Thanks to 
Lemma 3.6 again, this control can be directly transferred from the convexified measure v to the 
original measure fx. Taking these ingredients from [8] into account, the proof of Lemma 3.6 of [7] 
goes through for the non-convex case as well. 

Finally, we have a very weak control that holds for all points (see (1.7) in [8]): for any C > 
there are positive constants Co , ci and C2 such that 

P^(|Afc-7fcl >C^)<C^oexp(-ciiV=^). (4.3) 

Given C, we can choose the positive constants Co,ci,C2 to be the same in (4.1), (4.2) and (4.3). 

The set TZl,k in (3.7) was exactly defined as the set of events that these three rigidity estimates 
hold. From (4.1), (4.2) and (4.3) we have 

F^{nL^K{'^,a)) > l-Cocxp(-ciiV^=) (4.4) 

for any > 0, a > with some positive constants Co,ci,C2 that depend on i> and a. 

Lemma 4.1 Let L and K satisfy (3.1) and 6 is the exponent appearing in (3.1). Then for any 
small ^ and a there exists a set TZ* = 7^^_^_^(^2^/2, a/2) C TZl,k{^'^S/2, a/2) such that 

P''(7^*) > l-Coexp(- iciiV"^) (4.5) 

with the constants Co,Ci,C2 from (4.1). Moreover, for any y G TZ* we have 

\E^-Xk -ak\< CN-^K^, k e Il,k, (4.6) 
where au was defined in (3.4). 
Proof. For any u > define 

T^*L,Kj^,c,):^{yeTZ(iy,a) : P^- - 7fc| > iV-i+'') < cxp ( - ici7V^=) , e Il,k} (4.7) 

with the i/-dependent constants ci,C2 > from (4.4). Note that TZ* , unlike TZ, depends on the 
underlying measure fi through the family of its conditional measures fiy. Applying (4.4) for v = 
and setting TZ = TZL,K{C^S/2,a/2), TZ* = TZ*^ j^^^{C^S/2,a/2), we have 

PM(7^*) > i_Coexp(- iciiV^^) 
with some Cq, ci, C2. Now if y e 7?.*, then 

lE^-Xfc -7fe| < Coe-"i^''/3 + CAr-iiC«', fc e Il.k- (4.8) 
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In order to prove (4.6), it remains to show that \ak — 7fc| is bounded by CN ^K^ for any k G Il,k- 
To see this, we can use that g £ away from the edge, thus 

g{x) = g{y) + 0{x - y) 

(recah that y is the midpoint of J). By Taylor expansion we have 

k~{L~K~l)^N g = N Q+0{Ni^''^)=N\-ik-yL-K-i\Q{y)+0{N\j\'+N^'/^), 

•llL-K-1 JyL-K-1 

i.e. 

7fc = VL-K^i + ^ ~ ^^^^^^ ^ ^ + 0{N-^K^). (4.9) 

Here we used that J = Jy satisfies (3.22) according to Lemma 3.5, since y e T^l j<-(^^(5/2, a/2) C 
'^L,i<'(C'^/2, a). Comparing (4.9) with the definition of a^, (3.4), using (3.22) and the fact that 
y~ VL-K-i = ^\J\, we have 

\ak-lk\<CN-^Ki. (4.10) 
Together with (4.8) this impUes (4.6) and this completes the proof of Lemma 4.1. q 

4.2 Completing the proof of Theorem 2.3 

Wc first notice that it is sufficient to prove Theorem 2.3 for the special case m ~ N/2, i.e. when 
the local statistics for the Gaussian measure is considered at the central point of the spectrum. 
Indeed, once Theorem 2.3 is proved for any V, k and m = N/2, then with the choice V{x) = x^/2 
we can use it to establish that the local statistics for the Gaussian measure around any fixed index 
k in the bulk coincide with the local statistics in the middle. So from now on we assume m ~ N/2, 
but wc carry the notation m for simplicity. 

Given k and m = N/2 as in (2.16), we first choose L,L,K, satisfying (3.1) (maybe with a 
smaller a than given in Theorem 2.3), so that k ~ L + p, m = L + p hold for some \p\ < K/2. In 
particular 

\L- N/2\<K. (4.11) 

For brevity, we use fi — fiy and ]1 ~ j/^ in accordance with the notation of Theorem 3.1. 

We consider y e R*^ ^- ^{^"^5/2, a) and y e i?~ ^ ^{^'^5/2, a), where 6 is the exponent appearing 

in (3.1). We omit the arguments and recall that 

A^(^I,A>)>l-Cocxp(-iciiV^^), M(i?~^,^)>l-Coexp(-ici7V^^) (4.12) 
with some positive constants. 

Proposition 4.2 With the above choice of the parameters and for any y G i?2 ^{S.'^S/2, a) and 
y & R~ ,,JE'^5/2,a), we have 

E''yO{{NgY+j^){xL+p - xl+^+i), • • . {NgY+p){xL+p - ^L+p+n)) (4.13) 



E^W{{Ng1 ){x~,^^ - x~,^^^,), . . . {Ngf )ix~,_^^ - x~,^^^J) 



-L+p'^ i^+P L+p+\n---\ ^L+p' 

where e is from Theorem 3.1. 
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Theorem 2.3 follows immediately from (4.12) and this proposition. q 
The rest of this section is devoted to the proof of Proposition 4.2. 

Proof of Proposition 4. 2. We will apply Theorem 3.1, but first we have to bring the two 
measures onto the same configuration interval J to satisfy (3.8). This will be done in three steps. 
First, using the scale invariance of the Gaussian log-gas we rescale it so that the local density 
approximately matches with that of ■ This will guarantee that the two configuration intervals 
have almost the same length. In the second step we adjust the local Gaussian log-gas fly so that 
Jy has exactly the correct length. Finally, we shift the two intervals so that they coincide. This 
allows us to apply Theorem 3.1 to conclude the local statistics are identical. 

The local densities qv around ^L+p.v and qg around 7£^p q may considerably differ. So in the 
first step we rescale the Gaussian log-gas so that 

Qv{lL+p,v) = £'g(7Z+p,g)- (4-14) 

To do that, recall that we defined the Gaussian log-gas with the standard V{x) = external 
potential, but we could choose Vs{x) = with any fixed s > and consider the Gaussian 

log-gas 

^ cxp ( - NpnsW), H.(A) := - ^ y,(AO - _ ^ log I A, - A, I . 

i=l i<j 

This results in a rcscaling of the semicircle density go to gQ{x) := sgdsx) and 7i^G to q := 
s~Wi,G for any i, so gaili^G) gets rescaled to gai'^lc) = sgchi^c)- In particular, £'g(7l+p,g) 
is rescaled to s^dlz+p g"^^ ^^'^ thus choosing s appropriately, we can achieve that (4.14) holds 
(keeping the left hand side fixed). Set 

Os{x) -.^ 0{{Ng-^{'y^ Q)){x„i ~ X,n+l), ■ ■ ■ , 9Gil?n,G))i^ni -Xm+n)), m= L+p, 

and notice that C's(x) = 0(sx). This means that the local gap statistics &gOs is independent 
of the scaling parameter s, since the product {N g'^){xm — Xm+a) (notation defined in (2.15)) is 
unchanged under the scaling. So we can work with the rescaled Gaussian measure. For notational 
simplicity wc will not carry the s parameter further and we just assume that (4.14) holds with the 
original Gaussian V{x) — .t^/2. 

We have now achieved that the two densities at at some points of the configuration intervals 
coincide, but the lengths of these two intervals still slightly differ. In the second step we match 
them exactly. Since y e Rl^k{£,5 /2,a) and y £ Rj^ ^(^5/2, a), from (3.22) in Lemma 3.5 we see 
that 

I Jyl = IVL+K+i - VL-K-il = + 0{N-^K^) (4.15) 



Ngviv) 
JC 

NgGiy) 



I Jyl = IvL+K+i - VL-K-il = + OiN-'K^). (4.16) 



Since gy is C^, we have for any \j\ < K, 



\Qviy) - gviiL+j.v)\ <C\y - -fL+jy\ 

<C\y - vlM + C\jL+jy - IlM + 0{N-^K^) < CKN-\ 
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and similarly for gG{y). 

Using (4.15), (4.14) and that the densities are separated away from zero, we easily obtain that 



1^ satisfies s = Sy,y = 1 + 0{K-^+^). (4.17) 



Note that this s is different from the scaling parameter in the first step but it will play a similar 
role so we use the same notation. For each fixed y, y we can now scale the conditional Gaussian 
log-gas /Xy by a factor s, i.e. change y to sy, so that after rescaling | Jy| = l^syl- 
We will now show that this rescaling does not alter the gap statistics: 

Lemma 4.3 Suppose that s satisfies (4.17) and let fi — ^iq be the Gaussian log-gas. Then we have 

I [E''* - E^-'*]C'(x)| < CK-^+^ (4.18) 

with 

0{x) 0{{Ng'^){x^ - x^+i), . . . , {Ng^){xm - Xm+n)) 
for any L — K<m<L + K — n (note that the observable is not reseated). 
Proof. Define the Gaussian log-gas 

with T-L~ defined exactly as in (3.6) but Vy{x) is replaced with 

lf(x) = y,(x)-^^log|a;-y,|, Vs{x)^\fx\ I := {L - + K}. 

Then by scaling 

W'W{x) = E^sC'(x/s) = W^O{yi) + 0{\s - 1|), (4.19) 

where in the last step we used that the observable O is a smooth function with compact support. 
The error term is negligible by (4.17) and (3.1). 

In order to control [E^y — E^y]C'(x), it is sufficient to bound the relative entropy S{^~\^y). 
However, for any y G Rl,k wc have 

K > min 1 V , ^ > ^ (4.20) 

with a positive constant. Applying this for y, we see that satisfies the logarithmic Sobolev 
inequality (LSI) 

CK 

5(M|lMy) < —D{^,l\^,y), 

where 



is the relative entropy and the relative Dirichlct form of two probability measures. Therefore 
^(A^llMy) < j^^"^ E - ^^'(^01' = CK{s^ ifV.^^ CK^N-^s if 
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In the last step we used (4.11) which, by rigidity for the Gaussian log-gas, guarantees that \xi\ < 
CK/N with very high probability for any i E I. Together with (4.19) and (4.17) we obtain (4.18) 
and this proves Lemma 4.3. q 

Summarizing, we can from now on assume that (4.14) holds and that y, y satisfy | Jy | — \Jy\. By 
a straightforward shift we can also assume that Jy = Jy so that the condition (3.8) of Theorem 3.1 
is satisfied. The condition (3.9) has already been proved to hold in Lemma 3.5. Condition (3.10) 
follows from the definition of the sets TZ*t ,^ ,, and TZ- see Lemma 4.1. Thus all conditions 

of Theorem 3.1 are verified. Finally, we remark that the multiplicative factors 0i,+p ^Z+p 
(4.13) coincide by (4.14) and (2.15). Then Theorem 3.1 (with an observable O rescaled by the 
common factor fii^+p = Qj ) implies Proposition 4.2. q 



5 Gap universality for Wigner matrices: proof of Theo- 
rem 2.2 

In our recent results on the universality of Wigner matrices [21, 27, 28], we established the univer- 
sality for Gaussian divisible matrices by establishing the local ergodicity of the Dyson Brownian 
motion (DBM). By local ergodicity we meant an effective estimate on the time to equilibrium for 
local average of observables depending on the gap. In fact, we gave an almost optimal estimate 
on this time. Then we used the Green function comparison theorem to connect Gaussian divisible 
matrices to general Wigner matrices. The local ergodicity of DBM was done by studying the flow 
of the global Dirichlet form. The estimate on the global Dirichlet form in all these works was suffi- 
ciently strong so that it implied the "ergodicity for locally averaged observables" without having to 
go through the local equilibrium measures. In the earlier work [23], however, we used an approach 
common in the hydrodynamieal limits by studying the properties of local equilibrium measures. 
Since by Theorem 3.1 we now know the local equilibrium measures very well, we will now combine 
the virtue of both methods to prove Theorem 2.2. To explain the new method we will be using, we 
first recall the standard approach to the universality from [21, 27, 28] that consists of the following 
three steps: 

i) rigidity estimates on the precise location of the eigenvalues. 

ii) Dirichlet form estimates and local ergodicity of DBM. 

iii) Green function comparison theorem to remove the small Gaussian convolution. 

In order to prove the single gap universality, we will need to apply a similar strategy for the 
local equilibrium measure /Xy. However, apart from establishing rigidity for /iy, we will need to 
strengthen Step ii). The idea is to use Dirichlet form estimates as in the previous approach, but 
we then use these estimates to show that the "local structure" after the evolution of the DBM for 
a short time is characterized by the local equilibrium ^y in a strong sense, i.e. without averaging. 
Since Theorem 3.1 provides a single gap universality for the local equilibrium /Zy, this proves 
the single gap universality after a short time DBM evolution and thus obtain the strong form of 
the Step ii) without averaging the observables. Notice that the key input here is Theorem 3.1 
which contains an effective estimate on the time to equilibrium for each single gap. We will call 
this property the strong local ergodicity of DBM. In particular, our result shows that the local 
averaging taken in our previous works is not essential. 
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We now recall the rigidity estimate which asserts that the eigenvalues Ai, A2, . . . , Aat of a gen- 
eralized Wigner matrix follow the Wigner semicircle law qg{x) (2.6) in a very strong local sense. 
More precisely, Theorem 2.2. of [26] states that the eigenvalues are near their classical locations, 
{7j}^ii (2.7), in the sense that 

pJbj : |A, -7,1 > (logiV)^[min(j,iV-j + l)l"'^V-2/3l <Cexp[-c(log7V)«] (5.1) 



for any exponent C, satisfying 



Aq log log iV < C < 



log(lOA^) 
10 log log iV 



where the positive constants C, Aq, depend only on Ci„/, Cgup, 6*1, ^2, see (2.2), (2.4). In partic- 
ular, for any fixed a > and ;/ > 0, there are constants Co, Ci, C2 > such that for any N >\ and 

k e {aN, (1 - a)Nl, we have 

P (|Afc - 7fc| > N-^+'') < Co exp ( - ciiV^^) (5.2) 

and (5.1) also implies 

N 

" Ik f < N-^+^'' (5.3) 



E^(A, 



fc=i 

for any v > {). The constants Co,ci,C2 may be different from the ones in (4.1) but they play 
a similar role so we keep their notation. With a slight abuse of notation, we introduce the set 
T^L.K = T^L. ct) from (3.7) in the generalized Wigner setup as well, just 7/,. denote the classical 
localitions with respect to the semicircle law, see (2.7). In particular (5.1) implies that for any 

> 

¥{TZl,k{C, «)) > 1 - Co exp ( - ciN'-) (5.4) 

holds with some positive constants Co,ci,C2, analogously to (4.4). We remark that the rigidity 
bound (5.1) for the generalized Wigner matrices is optimal throughout the spectrum and it gives 
a stronger control than the estimate used in the intermediate regime in the second line of the 
definition (3.7). For the forthcoming argument the weaker estimates are sufficient, so for notational 
simplicity wc will not modify the definition of TZ. 

The Dyson Brownian motion (DBM) describes the evolution of the eigenvalues of a flow of 
Wigner matrices, H = Ht, if each matrix element hij evolves according to independent (up to 
symmetry restriction) Brownian motions. The dynamics of the matrix elements arc given by an 
Ornstein-Uhlenbeck (OU) process which leaves the standard Gaussian distribution invariant. In 
the Hermitian case, the OU process for the rescaled matrix elements Vij N^^^hij is given by the 
stochastic differential equation 

du,j = d/3ij - iftjdt, i,j^l,2,...N, (5.5) 

where /3ij, i < j, are independent complex Brownian motions with variance one and f3ii are real 
Brownian motions of the same variance. The real symmetric case is analogous, just f3ij are real 
Brownian motions. 
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Denote the distribution of the eigenvalues A = (Ai, A2, . . . , A^r) of Ht at time t by /t(A)/i(dA) 
where the Gaussian measure /i is given by (2.5). The density ft = ft.N satisfies the forward 
equation 

dtft=Cft, (5.6) 



where 

with /3 = 1 for the real symmetric case and /3 = 2 in the complex hermitian case. The initial data 
/o given by the original generalized Wigner matrix. 

Now we define a useful technical tool that was first introduced in [21]. For any t > denote 
by = an auxiliary potential defined by 



^ 1 
M/^(A):=^M//(A,), W][\):=-{\,-^,)\ (5.8) 

i=i 

i.e. it is a quadratic confinement on scale \/t for each eigenvalue near its classical location, where 
the parameter t > will be chosen later. 

Definition 5.1 We define the probability measure df-i'^ := Z^^e^^^'^ , where the total Hamilto- 
nian is given by 

'H^:=n + W^. (5.9) 

Here % is the Gaussian Hamiltonian given by (2.5) and Zt = Z^t is the partition function. The 
measure /j,"^ will be referred to as the relaxation measure with relaxation time t. 

Denote by Q the following quantity 

N 

0<t<l J ^ 
J = l 

Since is a generalized Wigner matrix for all i, (5.3) implies that 

Q < N-^+^'' (5.11) 

for any > if > A'o(j^) is large enough. 

Recall the definition of the Dirichlet form w.r.t. a probability measure ui 



N 



and the definition of the relative entropy of two probability measures guj and lo 

S{gu}\u}) := J gloggduj. 
Now we recall Theorem 2.5 from [291: 
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Theorem 5.2 For any t > and consider the local relaxation measure n'^ . Set ip :— and let 
9t '■= ft/''P- Suppose there is a constant m such that 

< CiV". (5.12) 

Then for any t > tN^ the entropy and the Dirichlet form satisfy the estimates: 

S{ftf^'\fin<CN^QT-\ D^\^t)<CN^QT-^, (5.13) 

where the constants depend on e' and m. 

Corollary 5.3 Fix a > and let t > N^". Under the assumptions of Theorem 5.2, for any 
t > tN^ the entropy and the Dirichlet form satisfy the estimates: 

D^iVIt) < CN^Qt^^. (5.14) 

Furthermore, if the initial data of the DBM, fo, is given by a generalized Wigner ensemble, then 

D^i^t) < CA^2"+2'^ (5.15) 

for any h' > 0. 

Proof. Since gt — ft/i^, we have 

N 

i=i •' 

N 

< 2^D''^{y/Yt) + 2N^QT~'^. 



i=l 



Thus (5.14) follows from Theorem 5.2. Finally, (5.15) follows from (5.14) and (5.11). g 

Define fy to be the conditional density of ffi given y w.r.t. /iy, i.e. it is defined by the 
relation /y/iy = {flj)y- For any y g TZl,k wc have the convexity bound (4.20). Thus we have the 
logarithmic Sobolev inequality 

^(/yMyl/^y) < C^Y. ^^Ty) (5.16) 



N 



and the bound 



To control the Dirichlet forms Di for most external configurations y, we need the following Lemma. 



21 



Lemma 5.4 Fix a > 0, > 0, and r > iV^°. Suppose the initial data fo of the DBM is given by 
a generalized Wigner ensemble. Then, with some small e' > 0, for any t > tN^ there exists a set 
Gl,k C TZl,k of good boundary conditions y with 

P^'''{gL,K)>l-CN~'\ (5.18) 

such that for any y G Gl,k we have 

E^^'''(v^)^C'iV'''™''' ft,y = {ft)y, I = Il,k, (5.19) 
iei 

and for any bounded observable O 

|[]E/t,yPy _E^^]0(x)| < C'i^l/2^2e' + a+i.-l/2^ (5 20) 

Furthermore, for any k Cz I we also have 

\Ef'■y^'yxk - 7fc| < CN-^+r (5.21) 

Proof. In this proof, we omit the subscript t, i.e. we use f ^ ft- By definition of the conditional 
measure and by (5.15), we have 

iGl iel 

By Markov inequality (5.19) holds for all y in a set §1^^ with f^^'{Ql,K) > 1 - CiV"^^'- Without 
loss of generality, by (5.4) we can assume that Q\ ^ C TZl.k- The estimate (5.20) now follows 
from (5.19) and (5.16). 

Similarly, the rigidity bound (5.2) with respect to //x can be translated to the measure /y/iy, 
i.e. there exists a set Q\ ^, with 

P^^(ai,K) >l-^?oexp(-ici7V^^), 
such that for any y g Cj£ and for any k £ I , we have 

ff-^'y(\xk - 7fc| > iV"'+") < exp ( - ^ciTV^^). 

In particular, we can conclude (5.21) for any y ^ Gj^ Setting Gl,k ■= Q\ k ^ k proved 
the lemma. q 

Lemma 5.5 Fix a > 0, v > 0, and r > N^°. Suppose the initial data /o of the DBM is given 
by a generalized Wigner ensemble. Then, with some small e' > 0, for any t > tN^ , k £ I and 
y G Ql,k, we have 

\Wyxk - E^'-^'^^Xfcl < A"iV-^/2+.^+»+2e; (-5 22) 

In particular, if the parameters chosen such that 

j^^-3./2+u+a+2e' < ^-^K^, and N-^ + '' < N'^ 

with some small ^ > 0, then 

\E>'yxk-ak\<CN~^K^, kel, (5.23) 
where au is defined in (3.4). In other words, the analogue o/(3.10) is satisfied. 
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Notice that if we apply (5.20) with the special choice 0(x) = Xk then the error estimate would 
be much worse than (5.22). We wish to emphasize that (5.23) is not an obvious fact although we 
know that it holds for y with high probability w.r.t. the equilibrium measure /x. The key point of 
(5.23) is that it holds for any y G Gl,k and thus with "high probability" w.r.t /t/i! 

Proof. Once again, we omit the subscript t. The estimate (5.23) is a simple consequence of 
(5.22), (5.21) and (4.10). To prove (5.22), we run the reversible dynamics 



(5.24) 



starting from initial data qo = fy, where the generator Cy is the unique reversible generator with 
the Dirichlct form D^^, i.e.. 



iei 



Recall that from the convexity bound (4.20), tk = K/N is the time to equilibrium of this dynamics. 
After differentiation and integration wc get, 



E^-Xfe - E-^-"-Xfc| = I ^ du^ Jidkqu)dfiy\ + 0(exp {-cK^')). 
From the Schwarz inequality with a free parameter R, we can bound the last line by 



1 

N 



du 



RidkVq^y + R-']d^iy + 0(exp i-cK' )). 



Dropping the trivial subexponential error term and using that the time integral of the Dirichlet 
form is bounded by the initial entropy, wc can bound the last line by 



RS{fy^ly\^iy) 



NR 



Using the logarithmic Sobolcv inequality for ^y and optimizing the parameter R, we can bound 
the last term by 



Xk 



Xk 



< tkR i^r^ ( Viy) + + 0(exp (-CA-' )) 



< 



iVl/2 



NR 

1/2 
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E^rny^) +0(exp(-cif^')). 



(5.25) 



Combining this bound with (5.19), wc obtain (5.22). q 

Wc now prove the following comparison for the local statistics of fi and ftfi^ where fi is the 
Gaussian /3-ensemble, (2.11), with quadratic V, and ft is the solution of (5.6) with initial data fo 
given by the original generalized Wigner matrix. 

Lemma 5.6 Fix 7i > 0, a > and t > . Then for sufficient small a there exist positive e and 
e' such that for any t > tN^ , for any n and for any compactly supported smooth observable O we 
have 



[E^'f - E"] 0(iV(x, - x,+i), N{xj - Xj+2), . . . , N{xj 
for any j G laN, (1 — a)A^] and for any sufficiently large N . 



)) 



< CN' 



(5.26) 
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Proof. We will apply Lemma 5.4, and we clioose L = j. Since K < N^^^, the right hand side 
of (5.20) is smaller than N^^ . Then we have 



< CN- 



(5.27) 



for all y g Gl.k with the probability of Qlj< satisfying (5.18). Choose any y S Tl* defined in 
Lemma 4.1. We now apply Theorem 3.1 with both /iy and /iy given by local Gaussian /3-ensemble. 
Thus the estimate (3.10) is guaranteed by (4.6) and (5.23). Since y,y G 7^ = TZL.K{£,^S/2,a) and 
Lemma 3.5 guarantees (3.9), we can apply Theorem 3.1 so that 



[E^y - E^y] 0(iV(xj - xj+i), N{xj - xj+2), . . . , ^(a;^ 



)) 



< CN- 



(5.28) 



for ah y e Ql.k and y G U*. Since P^(7e*) > 1 - N-", see (4.5), we have thus proved that 



[E^- - E^]0(Ar(a;, - x,+i),N{xj - x.+a), . . . ,iV(x, - x,+„)) 



(5.29) 



for all y G Gl.k- From (5.27), (5.29) and the probability estimate (5.18) for Gl,k, with possibly 
reducing e so that e < e', we obtain that 



[E/tA. - E^]0{Nix, - x,+i), N{x, - x-,+2), . . . , N{x, ~ z,+„)) 



< CN-". (5.30) 



This proves Lemma 5.6. q 

Recall that Ht is the generalized Wigner matrix whose matrix elements evolve by independent 
OU processes. Thus in Lemma 5.6 we have proved that the local statistics of Ht , for t > N^'^'^^'^ , 
is the same as the corresponding Gaussian one for any initial generalized matrix Hq. Finally, we 
need to approximate the generalized Wigner ensembles by Gaussian divisible ones. The idea of 
approximation first appeared in [19] via a "reverse heat fiow" argument and was also used in [46] 
via a four moment theorem. We will follow the Green function comparison theorem of [27, 28] and 
in particular, the result in [37] since these results were formulated and proved for the generalized 
Wigner matrices. 

Theorem 1.10 from [37] implies that if the first four moments of two generalized Wigner en- 
sembles, H"^ and H"^ , are the same and a certain level repulsion estimate holds for one of the 
ensembles, say, H^, then we have 



lim 



-K^piNix, - Xj + i),N{Xj - Xj+2),. . . ,N{Xj - Xj+n)) = 0. 



(5.31) 



In [37] it was assumed that one of the ensemble satisfies a level repulsion estimate. This estimate, 
however, was needed only for eigenvector universality; following the argument in [37] one can check 
that no level repulsion was actually used for the eigenvalue universality. Roughly speaking, the 
proof of [37] was based on the Green function comparison theorem [27] which compares eigenvalues 
at a fixed energy. In order to convert fixed energy to a fixed eigenvalue index, one needs to know 
that the total number of eigenvalues up to a fixed energy is the same for the two ensembles. The 
total number of eigenvalues up to a fixed energy E can be expressed in terms of integration of 
imaginary part of the trace of Green functions, i.e.. 



dy Im Tr 



1 

H - {y + h]) 
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with an rj slightly smaller than Thus the basic idea of the Green function comparison 

theorem can be employed and this leads to (5.31). As in the Green function comparison theorem, 
the matching of the fourth moments can be relaxed to that the differences of the fourth moments 
are of order N~'^ for any c > 0. 

We now apply (5.31) with the choice being the generalized Wigner ensemble for which we 
wish to prove the universality and H"^ = Ht with t = N~'^ for some small c'. Following [28], we 
construct an auxiliary Wigner matrix Hq such that the first three moments of Ht and the original 
matrix are identical while the differences of the fourth moments are less than N^^ for some 
small c" > depending on c' (see Lemma 3.4 of [28]). This completes the proof of (2.8) showing 
that the local gap statistics with the same gap-label {j) is identical for the Wigner matrix and 
the Gaussian case. The proof of (2.9) follows now directly from Theorem 2.3 that, in particular, 
compares the local gap statistics for different gap labels (fc and m) in the Gaussian case. This 
completes the proof Theorem 2.2. r-i 



6 Rigidity and level repulsion of local measures 
6.1 Rigidity of /iyi proof of Theorem 3.2 

We will prove Theorem 3.2 using a method similar to the proof of Theorem 3.1 in [7]. Theorem 
3.1 of [7] was proved by a quite complicated argument involving induction on scales and the loop 
equation. The loop equation, however, requires analyticity of the potential and it cannot be applied 
to prove Theorem 3.2 for a local measure whose potential Vy is not analytic. We note, however, 
that in [7] the loop equation was used only to estimate the expected locations of the particles. Now 
this estimate is given as a condition by (3.12) and thus we can adapt the proof in [7] to the current 
setting. For later application, however, we will need a stronger form of the rigidity bound, namely 
we will establish that the tail of the gap distribution has a Gaussian decay. This stronger statement 
requires some modifications to the argument from [7] which therefore wc partially repeat here. We 
now introduce the notations needed to prove Theorem 3.2. 

Let be a continuously diffcrcntiable nonnegative function with 6* = on [—1, 1] and 0" > 1 
for |a;| > 1. We can take for example 6{x) ~ {x ~ l)^l2;>i + (a; + l)^l2;<-i in the following. 

For any m S {aN, {l~a)Nj and any integer 1 < M < aN, we denote /('"■*^) = [m-M, m+Mj 
and M = = 2M + 1. Let /y := ^/3. For any k, M with \k - L\ < K - M, define 

Let 

where Zy^rfi is a normalization constant. Choose an increasing sequence of integers, Mi < M2 < 
...< Ma such that Mi = K^, Ma = CK^'^"^ with a large constant C, and M^/M-^_i ~ K"^ 
(meaning that cK^ < M^/M^-i < CK^). We can choose the sequence such that A < Cr]~^. We 

set uj~f :~ Lu^'^^''^ and we study the rigidity properties of the measures w^i, cja-i, • ■ • in this 
order. Note that fiy = uja since y e TZl.k = T^L,K{^S/2,a) guarantees that \xi — Xj\ < \ Jy\ < 
CK/N, see (3.22), thus for M = Ma = CK^''^'^ the argument of 6 in (6.1) is smaller than 1, so 
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= in this case. We also introduce the notation 



k+M 



2M + 1 ^ 



j=k-M 



Definition 6.1 We say that /.jy has exponential rigidity on scale £ if there are constants C',c, such 
that the following bound holds 



P^^dxfc - ak\ >£ + uK^N~^) < Ce- 



u> 0, 



(6.2) 



for any k Cz I . 

First we prove that /iy has exponential rigidity on scale MaN^^. Starting from 7 = A, by 
the Herbst bound and the logarithmic Sobolev inequality for /iy with LSI constant of order K/N 
(5.16), we have for any fc € {L - K + Ma, L + K - MaI that 



> 



< ^-c{N/K)Nb^ 



6 > 0, 



> 



< Ce" 



(6.3) 
(6.4) 



Using the estimate (5.23) we have that 



Thus we obtain 



[Ma 



N I 



\ < Ce- 



(6.5) 



Since x 



< a^fc < ^fc+M ^^'^ ^^'^ Q^fc's are regular with spacing of order we get 

[M] _ \M] 



^ [A/] [M] 
Xk- oik< x'^^\i - a' 

and we also have a similar lower bound. Thus 



k-M - -''k+M "fe+M + i^J'iJV 



'( 



Xk - OLk 



> CMaN-^ + 



N I 



\ < Ce" 



(6.6) 



foranyfce lL-/l+2il/^, L+ii:-2M^], where we used that M^A > A'«. Iffce {L-K, L-K+2MAj, 
then 

Xk - ak < XL-K+2MA - OtL-K+2MA + Cil/^TV""^ 

and 

Xk-ak> Vl-k-1 -ak> -CMaN^^. 

Thus we have the estimate 

\xk - ttfcl < \xl-K+2Ma - C(L-K+2Ma \ + CMaN~^. 

Since (6.6) holds for the difference xl-k+2Ma ^ c(l-k+2AIa, 'we have that it holds for Xk ~ otk as 
well (with at most an adjustment of C) for any fcs \L — K,L — K + 2MaI- Similar argument 
holds for fc G {L + K- 2Ma, L + Kj. Thus we proved (6.6) for ah G IL-K^L + Kj, i.e. we 
showed exponential rigidity on scale MaN~^. 



Now we use an induction on scales and we show that if 
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(i) for any fee |L - A' + 71/-^, L + K - M^} we have 

P^- (l^L^^"' - a^k'"^ I > uK^N-^^ < Ce-™' , u > 0; (6.7) 

(ii) exponential rigidity holds on some scale M^N~^, 

P^^^lxfe - afel > CAf^A^^i +uA'«iV-i) < Ce'™', fee/, u > 0; (6.8) 

(iii) we have the entropy bound 

5(/iy|w-,) < Ce-"^^'^"'", (6.9) 

then (i)-(iii) also hold with 7 replaced by 7 — 1 as long as My_i > A'^. The iteration can be 
started from 7 = A, since (6.7) and (6.8) were proven in (6.5) and in (6.6) (even with a better 
bound), and (6.9) is trivial for j = A since oja = fJ'y 

We first notice that on any scale M^, the bound (6.7) implies (6.8) by the same argument as 
we concluded (6.6) for any fc G / from (6.5). So we can focus on proving (6.7) and (6.9) on the 
scale 

To prove (6.9) on scale Af^_i, notice that (6.8), with the choice u = M^K^^, implies 

P^^dxfc - ttfcl > CM^N-^) < Ce""*^'^"'', kel. (6.10) 

Since 

unless -Xj\> CM^.iA^-^AT^'' > CM^N-'^K'^, we have that the scale CM^N''^ is by a factor 
AT'' smaller than the scale of Xi — xj built into the definition of (fj^^^^i-^)^ see (6.1). But for 
I, j e i(k,M.,-i) T^g^^^ -Xj\< \x, - ai\ + \xj -aj\+ CAf^_iiV-^ Thus (/jC^'^^t-i) = unless 
we are on the event described in (6.10) at least for one k. Moreover, |V0'^'^'^^^"^-' (x)| < N'-^ for 
any configuration x in J. Thus, following the argument in Lemma 3.15 of [7], via the logarithmic 
Sobolev inequality for iiy, we get 

^(My 1^7-1) < CA'7V-iE^-|V,^('^-'^^-'-i)|2 < CN^e-'^'"'''"' < Ce-^*^'-i^"'\ (6.11) 

Here we used that the prefactor N'~^ can be absorbed in the exponent by using that M^K~^^ — 
M'^_J<-^'^ > A'2«-5f' = K'' > iV^ Here we have used ^ = 877 and Af^_i > K^. We will not need 
it here, but we note that the same bound on the opposite relative entropy, 

is also correct. Thus (6.9) for 7 — 1 is proved. 

Now we focus on proving (6.7) on the scale M-y_i. Set 1 < M' < M < K and fix an index 
k G I such that \k — L\ < AT — M. We state the following slightly generalized version of Lemma 
3.14 of [7] 

Lemma 6.2 For any integers 1 < M' < M < K, k e {L ~ K + M,L + K - M} and k' G 

Ik - M + M', k + M - M'i, we have 




< Ce- 
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Compared with Lemma 3.14 of [7], we first note that TV^ in Lemma 3.14 [7] is changed to K^^ 
due to that 0('='*^)(x) in (6.1) is defined with a K^^ factor instead of A^"^. Furthermore, here we 
allowed the center at the scale M' to be different from k. The only condition is that the interval 
Ik' - M', k' + M'l C Ik- M,k + M}. The proof of this lemma is identical to that of Lemma 3.14 
of [7]. 

In particular, for any 7 = 2, 3, ... A and with M' = M-^-i and AI = < K'^M^^i and with 
any choice of k^ efL - K + M^,L + K - M^j, fc^_i e IL - K + M^-i, L + K - M-y-ij, so that 
- M^^i, fc^_i + A/..,_i| C Ik^ - M^, k^ + Af^l, we get 



— 1 



[A/-, 



> 



N 



< Ce 



(6.12) 



The entropy bound (6.9) and the boundedness of imply that 



where M^K''^'^ > if^^-s,, > ^r, ^ ^^^y combine it with (3.12) to have 

|E"^.Tfc -ak\< CK^/N. 

The measure in (6.12) can also be changed to fiy at the expense of an entropy term S{fiy\ujy). 
Using (6.9), we thus have 



[M., 



[M. 



-,] ( 



[M-,- 



N 



< Ce- 



Ce 



Combining it with (6.7) and recalling ^ = 3?7, we get 



a I. 



\M., 



< Ce- 



Ce 



(6.13) 



(6.14) 



This gives (6.7) on scale M^_i if u < cM^K ^''/^ with a small constant c. Suppose now that 
u > cM-fK^^^/"^ , which, in particular, means that u > cK-''!'^. Then, by (6.8), we have 



< 



[A/-,_ 



> CM^N-^ + (1 - CK-'^^^)u^^ 



< ^P^'y (\xk - afcl > CM^N-^ + (1 - CR-i/^y 



N J 



This proves (6.7) for 7 — 1. Note that the constants slightly deterioriate at each iteration step, 
but the number of iterations is finite (of order 1/?/ = 3/^), so eventually the constants C, c in 
(3.13) may depend on ^. In fact, since the detcrioriation is minor, one can also prove (3.13) with 
^-independent constants, but for simplicity of the presentation we did not follow the change of 
these constants at each step. 
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After completing the iteration, ft-om (6.8) for 7=1, Mi = K^, we have 

P^^ (|.Tfc - afcl > CK^N-^ + uK^N-^^ < Ce"™', k e I; 

This concludes (3.13) for m > 1. Finally, (3.13) is trivial for u < 1 if the constant C is sufRciently 
large. This completes the proof of Theorem 3.2. q 

6.2 Level repulsion estimates of Hyi proof of Theorem 3.3 

We now prove the level repulsion estimate. Theorem 3.3, for the local log-gas fiy with good bound- 
ary conditions y. There are two key ideas in the following argument. We first recall the weak level 
repulsion estimate (4.11) in [7], which in the current notation asserts 

P'^^ixL-K - VL-K-l < S/N) < CNs 

for any s > 0, and similar estimates may be deduced for internal gaps. Compared with (3.14), this 
estimate does not contain any /? exponent, moreover, in order to obtain (3.17), the N factor has 
to be reduced to (neglecting the irrelevant log factor). Our first idea is to run this proof for 
a local measure with only particles to reduce the N factor to K^. The second idea involves 
introducing some auxiliary measures to catch some of the /3 related factors. We first introduce 
these two auxiliary measures which are slightly modified versions of the local equilibrium measures: 

IJ-O ^J■yfl = Z(3{xl-K - VL-K-l)''^ i^y] ^ll '.^ ^J.yA = ZiW~'^ ^ly, (6.15) 
W = {XL-K - yL-K-l){xL-K+l - VL-K-l), (6.16) 

where Zq, Zi are chosen for normalization. In other words, we drop the term {xl-k — Vl-k-i)^ 
from the measure ^y in /ig and we drop in ni. To estimate the upper gap, yL+K+i — xl+k, 
similar results will be needed when we drop the term (j/l+a'+i — xl+k)^ and the analogous version 
of W, but we will not state them explicitly. We first prove the following results which arc weaker 
than Theorem 3.3. 

Lemma 6.3 Let L and K satisfy (3.1) and consider the local equilibrium measure /iy defined in 
(3.5). 

i) Let he any fixed positive constants and let y G TZl,k{£,5 /'2-,ot). Then for any s > we 
have 

F^y [xL-K - yL-K-i < s/N] < C {Ks log Nf+^ , (6.17) 

and 

F^'-lxL-K+i - yL-K-1 < s/N] <CiKs\ogNf^+' . (6.18) 

ii) Let y be arbitrary with the only condition that \yi\ < C for all i. Then for any s > we 
have the weaker estimate 

P^^ [xL-K - VL-K-i < s/N] < [tTJ) ' (^-1^) 

/ (75^- \ 2/3+1 

P^-^ [xL-K+1 - VL-K-l < s/N] < [jjjj , .7 = 0, 1. (6.20) 
To prove Lemma 6.3, we first prove estimates even weaker than (6.17)-(6.20) for /iy and fiyj. 
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Lemma 6.4 Let L and K satisfy (3.1). 

i) Let ^, a be any fixed positive constants and let y € TZl.k = TIl.k{(.S/2, a), then we have for 
any s > 

^"H^L-K ~ VL-K-1 < s/N) < CKslogN, (6.21) 
Wy-{xL-K-yL-K-i<s/N)<CKs\ogN, j = 0,l- (6.22) 

ii) Let y he arbitrary with the only condition that \yi\ < C for all i. Then for any s > we 
have the weaker estimate 

CsK 

^^H^L-K - VL-K-i < s/N) < — — , (6.23) 

l^^y I 

P'^-HxL-K-yL^K-i^s/N) <^^, J = 0,1. (6.24) 

Proof. We will prove (6.21), the same proof with only change of notations works for (6.22) 
case as well. We will comment on this at the end of the proof. 

For notational simplicity, wc first shift the coordinates by S such that in the new coordinates 
y = 0, i.e. yL-K-i = —yL+K+i and J is symmetric to the origin. With the notation a := —yL-K-i 
and / = fL — K,L + A' | . wc first estimate the following quantity, for any < (/? < c (with a small 
constant) 

pa — acp 



Z^:=l ... dx n - ^jfe-""^ ^^^'^""^^ 

J J -a+aif i jgj- 

= (1 - ^)^W(i^-i)/2 J ...J\lw Hiw, - w,)^e 



where we set 



By definition. 



(1 - (ys) ^XL+j , = J]^ dxL+j dw = J]^ Aw^. (6.25) 

\3\<K \3\<K 



-Q Jl (y,-(l-^)«;,.)^- 

k<L~K-\ fe>L+X+l 

(6.26) 

For the smooth potential V ^ we have 

V{S+(\-^)wi))-Y{S + w,)\<C\^w,\ < ^ (6.27) 

with a constant depending on V, where wc have used \wj\ < a < CK/N which follows from 
I Jyl < CK/N due to y e 7^L,/^, see (3.22). 

Using (1 — ip)wj — yk > (1 — ^ yk) for L — 2K <k<L — K— 1 and the identity 



(1 - (p)wj -yk = {wj - yk) 



Wj - yk 



30 



for any k, we have 



k<L~K~l 



k<L-K-l 



n V- 



(fiWj 



(6.28) 



and a similar estimate holds for k > L + K + 1. After multiplying these estimates for all j 
1, 2, . . . , A', we thus have the bound 



> 



mill 

\w\<a 



n 



ipw 



, k<L-2K 



w -yk 



n [> 



k>L+2K 



Vk~w 



K 



(6.29) 



Recall that y G TZl.k, i-e. we have the rigidity bound for y with accuracy N~^K^ <C K/N ^ a, 
see (3.7), i.e. y^'s are regularly spaced on scale a or larger. Combining this with |w| < a < CK/N, 
we have 



E 



k<L-2K 



Hence 



n [i- 



w - yk 



(fW 



< C(pK log N. 



w - yk 



> 1 - CipK log N, 



(6.30) 



(6.31) 



k<L-2K 

and similar bounds hold for the k > L + 2K factors. Thus for any (f < c we get 

^ > 1 - C(/3K^ + A'2 \ogN)(p > 1 - CA'V(logiV). 

Now we choose ip := s/{aN) and recall a ^ K/N . Therefore the //y-probability of xl+i — yL> 
aip = s/N can be estimated by 

P^-(xL-K - yL-K-1 > s/N) > ^ > 1 - CKs{\ogN). 

for all sA'log sufficiently small. If sA'log is large, then (6.21) is automatically satisfied. This 
proves (6.21). 

In order to prove (6.23), we now drop the assumption y G TZl.k and replace it with \yi\ < C. 
Instead of (6.28), we now have 



n {{i-p)w,-ykf>{i-^Y'' n 

k<L-K~l k<L-K-l 

and a similar estimate holds for fc > A + A' + 1. We thus have the bound 

K 

e-^^^f{l - > 1 - CipNK. 



(6.32) 



Pf'-ixL-K - yL-K-1 > s/N) >^> 

^0 



(6.33) 



With the choice ip := s/{\Jy\N) this proves (6.23). 

The proof of (6.22) and (6.24) for fiy Q is very similar, just the k = L — K — 1 factor is missing 
from (6.26) in case of j = —AT. For fiy^i, two factors are missing. These modifi-cations do not alter 
the basic estimates. This concludes the proof of Lemma 6.4. q 
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Proof of Lemma 6.3. Recalling the definition of fiQ and setting X := xl-k — Vl-k-i for brevity, 
we have 

From (6.22) we have 

E'^«[1(X < s/N)X^] < C{s/NfKs\ogN 
and with the choice s = cK^^iXogN)^"^ in (6.22) we also have 



JMo X > > 1/2 



with some positive constant c. This implies that 



2 \NK\ogN^ 

We have thus proved that 

P^>'[X < s/N]< C{s/Nf K slog N {NK log Nf = C {Ks log Nf^'^ , (6.35) 

i.e. we obtained (6.17). 

For the proof of (6.18), we similarly use 

P^'-IXL-K+I - VL-K-l < S/N] = W^ln^] ■ (^-^^^ 

From (6.22) we have 

E^-'[l{xL-K+i-yL-K-i < s/N)W^] < {s/N)^'^F'''[xL-K-yL-K-i < s/N] < Cis/NfEslogN. 
By the same inequality and with the choice s — cA'^^(log A^)^^, we have 

(w> I > 1/2 

V - {NKlogNY ) - ' 

with some positive constant c. This implies that 



2 \ {NK log Nf 
We have thus proved that 

Wy[xL-K+i - VL-K-i < s/N] < C{s/Nff^KslogN {{NKlogNff = C {KslogNf^^^ , 

(6.37) 

which proves (6.18). Finally, (6.19) and (6.20) can be proved using (6.23) and (6.24). This 
completes the proof of Lemma 6.3. q 

Proof of Theorem 3.3. For a given i, define the index set 

I := [max(i -K^,L-K - 1), min(i + A'^, L + K + 1)] 
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to be the indices in a neighborhood of i. We further condition the measure fjLy on the points 

Zj Xj j e I" := Il.k \ I 

and we let ^y_z denote the conditional measure on the remaining x variables {xj : j € I}. Setting 
L' ^ i. K' = K^, from the rigidity estimate (3.16) we have (y, z) € TZ = TZL',K'{£,^S/2,a) with a 
very high probability w.r.t. fiy. We will now apply (6.17) to the measure /iy^z with a new 5' = 
and K' = . This ensures that the condition < K' is satisfied and by the remark after (3.1), 
the change of d affects only the threshold A^o- We obtain 

F'^y-^ixi - .T,+i < s/N] < C {K^slogNf'^'^ (6.38) 

with a high probability in z w.r.t. ^y. The subexponential lower bound on s, assumed in part ii) 
of Theorem 3.3, allows us to include the probability of the complement of TZ in the estimate, we 
thus have proved (3.17). Similar argument but with (6.17) replaced by (6.18) yields (3.18). 

To prove the weaker bounds (3.14), (3.15) for any s > 0, we may assume that L — K < i < L; 
i > L is treated similarly. Since y € 'R-l,k, we have | Jy| > cK/N. We consider two cases, either 
Xi — UL-K-i < c'K/N or Xi — ul-k-i > c' K/N with c' < c/2. In the first case, we condition on 
XL-K, ■ ■ ■ ,Xi and we apply (6.23) to the measure vi = fJ-y,XL_K.---Xi- The configuration interval of 
this measure has length at least cK/(2N), so we have 

C'Ks 

P"^ (x.+i -x.< s/N) < ^j^jj^^ < CNs. (6.39) 

In the second case, Xi—yL~K~i > c' A'/A^, we condition on x^+i, 0:^+2, ... XL+i<-. The corresponding 
measure, denoted by V2 = fJ'y,xi^i,...XL+KJ has a configuration interval of length at least c'K/N. 
We can now have the estimate (6.39) for V2. Putting these two estimates together, we have proved 
(3.14). Finally (3.15) can be proved in a similar way. This completes the proof of Theorem 3.3. 

□ 



6.3 Comparison of the local statistics of two local measures; proof of 
Theorem 3.1 

In this section, we start to compare gap distributions of two local log-gases on the same configu- 
ration interval but with different external potential and boundary conditions. We will express the 
differences of gap distributions between two measures in terms of random walks in time dependent 
random environments. From now on, we use microscopic coordinates and we relabel the indices so 
that the coordinates of Xj are j G I = {—K, ...,0,1,... K}, i.e. we set L = L = in the earlier 
notation. This will have the effect that the labelling of the external points y will not run from 1 to 
N, but from some L_ < to i+ > with L+ — L_ = N. The important input is that the index 
set / of the internal points is macroscopically separated away from the edges, i.e. \L±\ > aN. 

The local equilibrium measures and their Hamiltonians will be denoted by the same symbols, 
fly and "Hy, as before, but with a slight abuse of notations we redefine them now to the microscopic 
scaling. Hence we have two measures fiy = e~^^^ /Zy and fiy = e~^^^ /Zy, defined on the same 
configuration interval J = Jy = Jy with center y, which, for simplicity, we assumed y = 0. The 
local density at the center is p(0) > 0. The Hamiltonian is given by 

i<3 



33 



V^{x) := W(x/iV)-2^1og| 



yj\ 



(6.40) 



and Hy is defined in a similar way with V in (6.40) replaced with another external potential V. 
Recall also the assumption that V",V" > —C (2.10). We will need the rescaled version of the 
bomids (3.22), (3.23) and (3.24), i.e. 



Uv 



/c 



d+{x) 



' d^{x) 
inf V" c 



O 



N 



d{x) 



d{x) 

X e J, 



X £ J, 



where 



(6.41) 
(6.42) 
(6.43) 

(6.44) 



d[x) := min{|x - |a; - yK+i\} 

is the distance to the boundary and we redefined d±{x) as 

d_{x) := d{x) + g{0)K^, d+{x) :== max{|x - y-K-i\, \x - Vk+iW + giO)K^. 

The rescaled version of Lemma 3.5 states that (6.41), (6.42) and (6.43) hold for anyy € TZl.k{£,S/'^, 
where the set TZl.k, originally defined in (3.7), is expressed in microscopic coordinates. 
Wc also rewrite (3.10) in the microscopic coordinate as 



\E''^Xj-aj\ + \E''fi 



< CK^, 



where 



J 



' /C + f 

is the rescaled version of the definition given in (3.4), but we keep the same notation. 
The Dirichlet form is also redefined; in microscopic coordinates it is now given by 



iei iei 

Due to the rescaling, the LSI from (5.16) now takes the form, for y e TZl.k, 

S{g^ly\^ly)<CKD^y{^). 
Define the interpolating measures 

< ~ = re [0,1], 



(6.45) 
(6.46) 

(6.47) 

(6.48) 
(6.49) 



so that cUy ~ = ^y and - = jiy [Z^ is a normalization constant). This is again a local log-gas 
with Hamiltonian 



(6.50) 
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and external potential 

Vl^{x): = {l~r)Vy{x)+rVy{x) 
Vy{x): = NV{x/N)^2Y,\og{x-v,), 

Vy{x) : = NV{xlN) - 2^1og(x -y,). 

The Dirichlet for w.r.t. the measure = - is defined similarly to (6.47). 
For any bounded smooth function O with compact support 

dr{0{Xp - Xp+i, ...Xp- Xp+n))uj^^ = l3(ho;0{Xp - Xp+i, ...Xp- Xp+n))u-^^, (6.51) 

where 

ho^ho{^)^Y.^Vy{xi)~Vy{x,)) (6.52) 

and {f]g)u '■= fg — (1E'^/)(IE'^.9) denotes the correlation. From now on, we will fix r. Our main 
result is the following estimate on the gap correlation function. 

Theorem 6.5 Consider two smooth potentials V,V with V",V" > — C and two boundary condi- 
tions, y,y G 7?.l=o.a:('?^'5/2, a), with some sufficiently small ^, such that J = Jy = Jy. Assume 
that (6.45) holds for both boundary conditions y,y. Then, in particular, the rescaled version of 
the rigidity bound (3.13) and the level repulsion bounds (3.17), (3.18) hold for both Hy and Jly by 
Theorem 3.2 and Theorem 3.3. 

Fix ^* > 0. Then there exist £ > and C > 0, depending on ^* , such that for any sufficiently 
small ^, for any < r < 1 and for \p\ < K^^^ we have 

|(/io; 0(xp - Xp+i, ...Xp- Xp+n))ui^\ < K^^K~^ (6.53) 

for any smooth function O : R" — > K with compact support provided that K > Ko{^,^* ,n) is large 
enough. 

Notice that this theorem is formulated in terms of K being the only large parameter; N 
disappeared. We also remark that the restriction |p| < K^^^ can be easily relaxed to \p\ < 
K — K^~^ with an additional argument conditioning on set {xi : i E I \ 1} to ensure that p is 
near the middle of the new index set /. We will not need this more general form in this paper. 

First we complete the proof of Theorem 3.1 assuming Theorem 6.5. 

Proof of Theorem 3.1. The family of measures ~, < r < 1, interpolate between /iy and /iy . 
So we can express the right hand side of (3.11), in the rescaled coordinates and with L = L = as 

[EA'y - Wi]0{Xp - Xp+i, ...Xp- Xp+n) 

Using (6.51) and (6.53) we obtain that this difference is bounded by K'~^^K^^. Choosing ^ suffi- 
ciently small so that K^^K^"^ < /v we obtain (3.11) (with e/2 instead of e). This completes 
the proof of Theorem 3.1. r-i 



< 



dr—K^yyOixp — a;„+i, 
dr 



Xp+ri)- (6.54) 
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In the rest of the paper we will prove Theorem 6.5. The main difficulty is due to the fact that 
the correlation function of the points, {xi\Xj)^, decays only logarithmically. In fact, for the GUE, 
Gustavsson proved that (Theorem 1.3 in [33]) 

N 

{xf,Xj}GUE log |j— -pp^: (6-55) 

and a similar formula is expected for w. Therefore, it is very difficult to prove Theorem 6.5 based 
on this slow logarithmic decay. We notice that, however, the correlation function of the type 

{9i{xi);g2{xj - Xj+i))^ (6.56) 

decays much faster in \i — j\ due to that the second factor g2{xj — Xj+i) depends only on the 
difference. Correlations of the form (5i(a;i — a;i+i);52(a;j —a;j+i))a; decay even faster. The fact that 
observables of differences of particles behave much nicer was a basic observation in our previous 
approach [21, 27, 28] of universality. 

The measure w = cj^ - is closely related to the measures and /iy. Our first task in Section 6.4 
is to show that both the rigidity and level repulsion estimates hold w.r.t. the measure uj. Then 
we will rewrite the correlation functions in terms of a random walk representation in Proposition 
7.1. The decay of correlation functions will be translated into a partial regularity property of the 
corresponding parabolic equation, whose proof will be the main content of Section 8. Section 7 
consists of various cutoff estimates to remove the singularity of the diffusion coefficients in the 
random walk representations. We emphasize that these cutoffs are critical at j3 = 1; we do not 
know if our argument can be extended to /3 < 1. 

6.4 Rigidity and level repulsion of the interpolating measure w^- 

In this section we establish rigidity and level repulsion results for the interpolating measure ~, 
similar to the ones established for /iy in Section 6 and stated in Theorems 3.2 and 3.3. 

Lemma 6.6 Let L and K satisfy (3.1) and y, y G TZl.k{S,^5 /2,a). With the notation = <^y y 
there exist constants C , 9^ , C2 and C3 such that the following estimates hold: 
i) [Rigidity bound] 

P"(|a;, -ai\> CK^'-^^) < Ce-^"', i G /. (6.57) 
a) [Weak form of level repulsion] For any s > we have 

r'{x,+i-x^<s)<C{Nsf^\ iefL- K -1,L + Kl s > 0, (6.58) 

r^{x,+2-x^<s)<C{Nsf^^\ ielL-K -l,L + K-ll s > 0, (6.59) 
Hi) [Strong form of level repulsion] With some small 9 > 0, for any s > exp {—K^) we have 

r^{x^+i-x^<s) <C{K'^''isf^\ i(zlL-K -1,L + Kl, (6.60) 

ir{x,+2-x^<s) <C{K^''isf^^\ ielL-K-l,L + K-lj, (6.61) 
iv) [Logarithmic Sobolev inequality] 

S{gu\u:) < CKD^i^). (6.62) 
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Note that in (6.57) we state only the weaker form of the rigidity bound, similar to (3.16). It is 
possible to prove the strong form of rigidity with Gaussian tail (3.13) for w, but we will not need 
it in this paper. 

The level repulsion bounds will mostly be used in the following estimates which trivially follow 
from (6.58)-(6.61): 

Corollary 6.7 Under the assumptions of Lemma 6.6, for any p < /3 + 1 we have 

-<CpA'^^«, i^lL-K-l,L + Kl (6.63) 

and for any p < 2/3 + I 

E" I — < CpK^'^, ielL-K -1,L + K -11 (6.64) 

□ 

The key to translate the rigidity estimate of the measures fiy and fiy to the measure uj — uj^ ~ 
is to show that the analogue of (6.45) holds for uj. 

Lemma 6.8 Let L and I'C satisfy (3.1) and y,y G TZl,k{S,5 /2,q). Consider the local equilibrium 
measure /iy defined in (3.6) and assume that (3.10) is satisfied. Let uj^ ~ he the measure defined in 
(6.49). Recall that denote the equidistant points in J, see (6.46). Then there exists a constant 
C, independent of such that 

E'^y.s \xj - ajl < CK^^. (6.65) 
Proof of Lemma 6. 8. We first prove the following estimate on the entropy. 

Lemma 6.9 Suppose pLi is a probability measure and ui — Z^^e^dfii for some function g and 
normalization Z . Then we can bound the entropy by 

S S{ijj\pLi) = E^g-\ogW^eS < E^^.g-E^i^. (6.66) 

Consider two probability measures d/i^ = Z~^e^^'dx, i = 1,2. Denote by g the function 

g = r{Hi - i?2), < r < 1 (6.67) 

and set lu = Z^^e^dfii as above. Then we can bound the entropy by 

min(S'(w|^i), 5(w|/^2)) < - E''^ 1 [Hi - H2). (6.68) 
Proof. The first inequality is a trivial consequence of the Jensen inequality 

S* = E"g - logE^ief < E'^^ - E^^^. 
The entropy inequality yields that 

E'^g <r\ogW^e^''' + rS. (6.69) 

By the definition of g, we have 

logE'^igf/'' = - logy e-sl'^A^l2 < W^g/r. 
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Using this inequality and (6.69) in (6.66), we have proved 



S < 



1-r 



{Hi - H2). 



(6.70) 



We can assume that r < 1/2 < 1 — r since otherwise we can switch the roles of Hi and H2. Hence 
(6.68) holds and this concludes the proof of Lemma 6.9. q 

We now apply this lemma with ^2 = My ^-^d /ii = Hy to prove that 

min[5(c.;~|^y),5(c.;-|My)] < A^^« 
To see this, by definition of g and the rigidity estimate (3.13), we have 



(6.71) 



W^g - W^g = 



As 



Vy{sai + (1 - s)xi) - V~(sa, + (1 - s)xi) {xi - Ui) 



[0,1] d{sai + (1 - s)xi) 



(6.72) 



In the first step we used that the leading term Vy{ai) — Vy{ai) in the Taylor expansion is determin- 
istic, so it vanishes after taking the difference of two expectations. In the last step we used that 
with a very high fii- or /^2-probability d{sai + (1 — s)xi) ^ d(ai) are equidistant up to an additive 
error if i is away from the boundary, i.e., + K'~^^ < i < K ~ K^^, see (3.13). For indices 
near the boundary, say —K < i < —K + K'~'^, wc used d{sai + (1 — 5)2;^) > cminjl, (i(a;_i4-)}. 
Noticing that d{x-K) = X-k — y-K-i, the level repulsion bound (3.17) (complemented with the 
weaker bound (6.17) that is valid for all s > 0) guarantees that the short distance singularity 
[d{x-K)]~~^ has an E^^'^ expectation that is bounded by CK^^. 

We now assume that (6.71) holds with the choice of S{ujy yl^y) for simplicity of notation. By 
the entropy inequality, we have 



E'^y-^\x, - a,\ < logE^^el^'-"'! +7^^^^. 
From the Gaussian tail of the rigidity estimate (3.13), we have 

logE^^el"^'-"'! < K'^^. 



(6.73) 



(6.74) 



Using this bound in (6.73) we have proved (6.65) and this concludes the proof of Lemma 6.8. q 
Proof of Lemma 6.6. Given (6.65), the proof of (6.57) follows the argument in the proof of 
Theorem 3.2, applying it to instead of ^. Once the rigidity bound (6.57) is proved, we can follow 
the proof of Theorem 3.3 to obtain all four level repulsion estimates, (6.58)-(6.61), analogously to 
the proofs of (3.14), (3.15), (3.17) and (3.18), respectively. The logiV factor can be incorporated 
into K^^^. 

Finally, to prove (6.62), let be the reversible generator given by the Dirichlet form 

^ / /^"/d-;,y = E / (6.75) 
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Thus for the Hamiltonian H = ~ of the measure uj = uj^ ~ (see (6.50)), we have 



by using (6.43) and d{x) < CK for good boundary conditions. Thus LSI takes the form 

S{gu\u) < CKD^i^). (6.77) 
This completes the proof of Lemma 6.6. q 



The dynamics given by the generator with respect to the interpolating measure w = w*^ 



can also be characterized by the following SDE 
dx, = dB,+P 



y.y 



where {B^k , B^k+i, ■ ■ ■ , Bk) is a family of independent standard Brownian motions. With a 
slight abuse of notations, when we talk about the process, we will use and to denote the 
probability and expectation w.r.t. this dynamics with initial data w, i.e., in equilibrium. This 
dynamical point of view gives rise to a representation for the correlation function (6.53) in terms 
random walks in random environment. 

Starting from Section 7 we will focus on proving Theorem 6.5. The proof is based on dynamical 
idea and it will be completed in Section 7.5. 

7 Local statistics of the interpolating measures: Proof of 
Theorem 6.5 

7.1 Random Walk Representation 

In this section we derive a random walk representation for the gap correlation function on the left 
hand side of (6.53). For brevity of the formulas, we consider only the single gap case, n = 1; the 
general case is a straightforward extension. The gap index p is fixed, —K < p < K — 1. 

Suppose h{t) = h{t,x.) is the solution of the equation dth = C^h where £ = is the imique 
reversible generator with Dirichlet form D'^{g) w.r.t. the measure lo = y- The initial condition 
is chosen to be /iq given in (6.52). Introduce the notation 

v(t,x) = Vx/i(t,x), i.e. wj(i,x) :=a,^/i(<,x), (7.1) 

then 

z;,(o,x) = a,^./jo(x) = (yy)'(x,)-(yy)'(x,), .?e/. (7.2) 

By integrating the time derivative of {h{t, x); 0{xp~Xp+i))u, and using the equation h{t^ x) satisfies, 
we have 

(/io(x);0(xp - Xp+i))^ = j da J 0'{xp- Xp+i)[vp{a,x) - Vp+i{cr,x)]duj{x). (7.3) 
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Now we derive a random walk representation for (7.3). Taking the gradient of the equation dth = 
Ch, a direct computation of the commutator [V,£] yields that 

atv(t, x) = C^^{t, x) - ^(x)v(t, x), (7.4) 



with initial condition vo(x) = v(0,x) = VhQ{x), see (7.2). Here A :^ B + W and B and W are 

the following x-dependent operators acting on vectors v G MJ^ with components u^, i e /: 

[S(x)v], = - ^ B,k{^){vk - Vj), B,fc(x) = ^ > (7.5) 

^ [Xj Xk) 

and 

[W(x)v], := W,-(x)z;, (7.6) 

with the function 



I |fe|>if+i 



1 -r 



l-r 



Here r e [0, 1] is a fixed parameter which we will omit from the notation of W. 

Recall that x(t) is the stochastic process (6.78) with generator C = and let Ex denote 
the expectation for this process with initial point x(0) = x. The expectation with respect to the 
process starting from equilibrium is E"[-] — /Ex[-]w(dx). The solution to (7.4) can be represented 
by the Feynman-Kac formula: 

v(a,x) =ExlSp(-^ ^(x(s))ds)vo(x(a)), (7.8) 

where 

E^p(-f A{x{s))ds) -.^ 1 - f ^(x(si))dsi+ / ^(x(si))^(x(s2))dsids2 + ... (7.9) 

is the time-ordered exponential. To prove that (7.8) indeed satisfies (7.4), we notice from the 
definition (7.9) that 

v(a,x) =Ex£^(-^ I(x(s))ds)vo(x(a)) (7.10) 
= ExVo(x(a)) Ex^(x(si))Ex(si)Sp(^ I(x(.s))ds) Vo(x(a))dsi. 
Using that the process is stationary in time, we have 

v((T, x) = ExVo(x(cr)) - / Ex^(x(si))v((7 - si, x(si))dsi 

"'0 

= ExVo(x(cr)) - / Ex^(x((T - si))v(si,x(ct - si))dsi 
Jo 

= e'^^vo(x) - /''[e('^-^i)^I(-)v(si, •)](x)dsi. 
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Differentiating this equation in a we obtain that v defined in (7.8) indeed satisfies (7.4). 

For any path x(-) starting from x fixed at time zero, x(0) = x, and for any vector wq, denote 
the vector ^ 

w(s;x(-),fT) :=lSp(-^ ^(x(6))d6)wo. (7.11) 

Then we have 

9sw(s;x(-),(t) = (9s 1- / ^(x(si))dsi + / ^(x(si))^(x(s2))dsids2 - • • ■ wo 

= ^(x(s))w(s;x(-),cr). 

This is the backward equation, starting from w(s = CT;x(-),cr) = Wg and evolving backwards in 
time from s = cr to s = 0, for a random walk given by the random environment induced by 
the path x(-). The equation in a slightly different setting already appeared in Proposition 2.2 of 
[16] (see also Proposition 3.1 in [32]), which was a probabilistic formulation of the idea of Helffer 
and Sjostrand [35] and Naddaf and Spencer [41]. For our purpose, we will convert the backward 
equation into the forward equation in order to be consistent with the convention in the partial 
differential equation literature [10]. Thus replacing s with cr — s in (7.11) we define 

w(s;x(-),cr) := w(cr - s;x(-),cr) = E^pf - /" ^(x(&))d6Wo. (7.12) 

J a — s 

We have then the evolution equation 

9sw(s;x(-),ct) = -^(x(cr - s))v^f(s;x(-),cr), v^f(0; x(-), cr) = Wo- (7.13) 

With these notations, and choosing wq = vo(x((t)) in the initial value problem (7.13), we have 
from (7.8) that 

v(a,x) =Exw(a;x(-),fT). (7.14) 

In other words, for any path starting from x and ending at x(cr), we solve (7.13), which is /C- 
dimensional system of ODE. Then taking expectation yields the solution v to (7.4). 

We now slightly rewrite the equation (7.13) into a more convenient form. Fix < it < 
CxKXogK and fix a path {x(s) : s e [0, a]}. Define the following operators on M'^ 

yl(s) ;=^(x(cr-s)), i3(s) := 6(x(a - s)), W(.s) W(x(cr - s)), (7.15) 

where W is a multipfication operator with the j-th diagonal Wj(s) = Wj{xj{a — s)) depending 
only the j-th component of the process x(s). Clearly A{s) = B{s) + yV(s). We also define the 
associated (time dependent) quadratic forms which we denote by the corresponding lower case 
letters, in particular 

b(.s)[u,v] = ^ii,[S(s)v]j = ^Y1 Bjk{s){uk - Uj){vk - Vj) 

iGl k.jel 

tT3(s)[u,v] = ^Ui[W{s)v]i =^W,Wi(s)Wj 
iel i 

and 

a(s)[u, v] - b(s)[u, v] + tt,(s)[u, v]. (7.16) 
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With this notation, (7.13) becomes 

a,w(s; x(-), a) = -^(s)w(s; x(-), a), w(0; x(-), a) = V/io(x(a)). (7.17) 

Combining (7.14) and (7.17) with (7.3), we have proved the foUowing proposition. Notice that 
in the above derivation we never used the exphcit form of the initial condition ho given in (6.52), 
so we can formulate the resuU for general ho- 

Proposition 7.1 For any smooth function /iq : J'^ — > R, for any p E I, —K < p < K — 1, we 

have 

{ho;0{xp - Xp+i))uj = J d<7 y 0'(xp - 2;p+i)Ex[wp(cr, x(-); cr) - Wp+i(cr, x(-); cr)]w(dx) (7.18) 

Here for any a > and for any path {x(s) G j'^ : s G [0, cr]} we let w be the solution of (7.17) 
with initial data w(0;x(-),(t) := V/io(x((t)). 

Now we apply this proposition to our case: 

Corollary 7.2 Let ho be given by (6.52) and assume that y,y G TZL=o,K{£,'^S/2,a). Then for any 
p € I, —K < p < K — 1, we have 

{ho;0{xp - Xp+i))^ (7.19) 
da / 0'{xp - Xp+i)E^[wp{a,x{-);a) - Wp+i{a,x{-);a)]uj{dx) + 0{K~^). 

Proof. Since y,y G TZl=o.k{£,'^^ f"^, ct), we easily obtain that Wj{s) > cK^^, which means that 
the L^-norm of the solution to (7.17) decays at a rate of order K. Therefore the integral in (7.18) 
from CiKlogK to infinity can be bounded by 

/•oo 

C / e-'='^/^E"|up(CT, ■)\da < CK^^'+^'>^K-^' < R-^ 

JciKlogK 

if the constant Ci is chosen sufficiently large. Here we used that 

E-\vp{a, .)| - E-\vp{0, .)| < E-[\{VyY{xp)\ + |(Fy)'(xp)|] (7.20) 



.,„,,, - 1 

< 



-\yj-xp\ Wj 



from (6.63) and using that y,y G TZl=o,k = T^l=o.,k{S,'^5 /2,a) are regular on scale K^^ < K^, so 
the summation is effectively restricted to terms. q 

The representation (7.18) and (7.19) expresses the correlation function in terms of the discrete 
spatial derivative of the solution to (7.17). To estimate Wp{a,x{-); a) — Wp+i{a,x{-);a) in (7.19), 
we will now study the Holder continuity of the solution w(s,x(-);(t) to (7.17) at time s = a and 
at the spatial point p. We will do it for each fixed path x(-), with the exception of a set of "bad" 
paths that will have a small probability. 

Notice that if all points Xi were approximately regularly spaced in the interval J, then the 
operator B has a kernel Bij ~ (i — j)~^, i.e. it is essentially a discrete version of the operator 
\p\ = V— A. Holder continuity will thus be the consequence of the De Giorgi-Nash-Moser bound 
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for the parabolic equation (7.17). However, we need to control the coefficients in this equation, 
which depend on the random walk x(-). 

For the De Giorgi-Nash-Moser theory we need both upper and lower bounds on the kernel Bij . 
The rigidity bound (6.57) guarantees a lower bound on Bij, up to a factor K~'-^'^^ > . The 
level repulsion estimate implies certain upper bounds on Bij. In the next section we define the 
good set of paths that satisfy both requirements. 



7.2 Sets of good paths 

From now on we assume the conditions of Theorem 6.5. In particular we are given some ^ > and 
we assume that the boundary conditions satisfy y,y G 'Rl=q.k = 'R-L=o,K{S,'^S/2,a) and (6.45) 
with this ^. We define the following "good sets" : 

g-.^l sup sup \xj{s) -aj\< K^'}, (7.21) 

0<s<Cii<'log_R- \j\<K 

where 

e' (C•2 + 1)C^ 

with C2 being the constant in (6.57) and aj is given by (6.46). For any Z E I and a G [0, Ci-R" log K] 
we also define the event 



Qa.z := I sup sup da ^ 

^ 0<s<a KKKK + S Jo "'^ 



m\Ti{M+Z,K} 

E ^ < K"], (7.22) 

where p > will be specified later and, by convention, Xi{a — a) ~ yi whenever |i| > K . Define a 
new good set Qa^z by 

Qa,Z = Qa^Z n Qa^-K H Qa,K, (7.23) 

i.e., this is the event that the gaps between particles near Z and near the boundary are not too 
small in an appropriate average sense. 

Lemma 7.3 There exists a positive constant 9, depending on ^' = (C2 + 1)^^, such that 

P"(g'=) < Ce"-^*. (7.24) 

Moreover, there is a constant C4, depending on the constant C2 in (6.57) and on C3 in (6.60), 
(6.61) such that for any ^ and p small enough, we have 

TP'^iQa.z) < CK^^^-P (7.25) 
for each fixed Z E I and fixed a £ [0,CiK log K] . 

Proof. From the stochastic differential equation of the dynamics (6.78) we have 



\x^it) - x^{s)\ <C|t-s| + 



/ \X]{a) - Xi{a)\^^Jyj - x^{a)\-^'^^^^'''*^ ^'^'^ 



(7.26) 
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Using (6.63) and that x(-) is invariant under w, we have the bound 



^ \x,(a 



(a) - Xi{a)\ 



3/2 



< C'j^3+C3?|^_g|3/2_ (7 27) 



This imphes for any fixed s < t < CiKlogK and for any R> that 

<CX3+C3C|i_s|3/2^-3/2^ 



(7.28) 



A similar bound holds for the second summation in (7.26); the summation over large j can be 
performed by using that y is regular, y G TZl=o,k- 

Set a parameter q < cR and choose a discrete set of increasing times {sk : ^ = 0,1, 2..., K^/q} 
such that 

- So < si < .S2 < ... < CiA'logA', and \sk - Sk+i\ < q. (7.29) 

From standard large deviation bounds on the Brownian motion increment Bi{t) — Bi{s) and from 
(7.26), we have the stochastic continuity estimate 

( sup \x,{s) - x,(t)\ >r\< Ke'^^'''^ + CK^q^'^R-^'^ 

for any fixed fc. Taking sup over fc, we have 

( sup \x,{s) - x,{t)\ > i? ) < K^q-^e-^'^'/'^ + CK%'/^R-^/^ (7.30) 

\a<s,t<CiK\ogK,\t-s\<q,\i\<K J 

for any positive q and R with q < cR. 

From the rigidity bound (6.57) wc know that for some ^3 > and for any fixed k we have. 



'[\xj{sk)-aj\>CK^-'i^^<Ce 



(7.31) 



Choosing R = K^' /2 and q = cxp (-AT''''/^), and using that CK^^^'^ < K^' /2 with the choice of 
we have 

P"(g=) < Ce-'^'^K^-^ + + ^^6^1/2^-3/2 < ^^^p (_^^e3/3)^ (7 

for sufficiently large K , and this proves (7.24) with 9 ~ 6*3/3. 

We will now prove (7.25). We will consider only the set Q% ^ ^"^^ only for Z = 0. The 
modification needed for the general case is only notational. We start the proof by noting that 



-^r'da^ f - 

1 + s' Jo M' ^j^^, \xi{a-a)- x^+l 



< C- 



1 + s 



(a-a)P 

M 



(7.33) 



da- y 



' M \xi{a - a) - Xj+i((T - a)\^ 
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holds for any s' £ [s/2, s] and M' £ [M/2, A/]. Hence it is enough to estimate the probabihty 

1 1 1 

/ da— y —. ^ rT^>KP] (7.34) 

ll + .s7o il/^^^Ja;,(cr-a)-x,+i(cr-a)|2 - J 

for fixed dyadic points {s,M) = {{2~pK^,2~'^K)} in space-time for each integer p,q < Clog A'. 
Since the cardinality of the set of these dyadic points is just C(log K)"^, it suffices to estimate (7.34) 
only for a fixed s, Af. 

The proof is different for [3 = 1 and /3 > 1. In the latter case, from (6.63) we see that the 
random variable in (7.34) has expectation CK'-^^^ . Thus the probability in (7.34) is bounded by 
CK^^^-B, so (7.25) holds in this case with d slightly larger than C3 + 1 to accommodate the 
C{\ogKf factor. 

The case /? = 1 is more complicated since the random variable in (7.34) has a divergent 
expectation. We will use stochastic continuity. The proof below would work for any /3 > 1 as well, 
but for simplicity we set /3 = 1. For any fixed time Sk from the set (7.29), for any i <E |— A' — 1, K\ 
and for any (5 > 0, we have 

P"{|x,(sfc) - x,+i(sfe)| < 3,5} < C{N5f 
from the weak form of level repulsion (6.58). Thus 

P"{|x.(sfc)-a;,+i(sfc)| < 3,5 : V ^ £ [-.A - 1, A'], Vfc = 0, 1, . . . A'Vg} (7.35) 
^ - 2^m(sfc)| < 3J} < CK^q-\N5f < NH\ 

i k 

Set 

Q.:=\ sup \xi{t) - Xi+i{t)\>5\. 

^ 0<t<C 1 K log K,\i\<K 

Combining (7.35) with the stochastic continuity (7.30) (with R = S), we get 

P" {n'') < CNH"^ + A'3q-ie-^*''/« + CA'^qi/2^-3/2_ 

Choosing 5 := exp (— A'^), q := 5^, we obtain 

P" (!]")< exp(-ci^«). (7.36) 

On the set J7 we can estimate (7.34) by Markov inequality. By the strong level repulsion (6.60), 
we have 



E 



M ^ 



1 + S Jo M .^3^^ |xi(cr - a) - a::j+i((7 - a)|2 



< c sup E" , < r 1 m < CK'^^i\\og5\ < CX(2C^+i)« 

0<t<s \x,,{t) - Xi+l{t)\^ Js £^ 

(notice that the condition 6 > exp(— AT^) set for (6.60) is satisfied for 9 = Therefore 

(7.34) < P"(fi'=) + CK'-^'^^+^'^^-p < CK'-^'^^+^'^^-p. 

Choosing C4 slightly larger than 2C3 + 1 (to accommodate the C(logA')2 factor from the dyadic 
argument), we have proved (7.25). q 
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7.2.1 Restriction to the set Q 

Now we show that the expectation (7.19) can be restricted to the good set Q with a small error. 
With a slight abuse of notations we use G also to denote the characteristic function of the set Q. 
We just estimate the complement as 

\0'{xp ~ Xp+i)\E,:^g''[wp{a,x{-); a) - Wp+i{a, a)]uj{dx) 

<C J E^g'[\wp{a,ic{-y,a)\ + \wp+i{cT,^{-);cT)Mdx). 

Since A>OasalCxlC matrix, the equation (7.17) is contraction in L^. Clearly ^ is a contraction 
in as well, hence it is a contraction in any i'', 1 < g < 2, by interpolation. By the Schwarz 
inequality and the ^''-contraction for some 1 < q < 2, we have 

E-g^\wp{aM-);<^)\ < [E'^G'Y^^'^'^E^iwpiaM-y^'^WV^' 

1/9 



< [F^g-Y^^" ''[E"^|u;,(0,x(.);a)r 



with some 64 > 0. Here we used (7.24) for the first factor. For the second factor we recall 
w(0, x(-); cr) = V/io(x(tT)), which can be estimated by (6.63) using the invariance of the measure 
ui (recall the definition of ho from (6.52)). 
Hence we have proved that 



|0'(.Tp -.Tp+i)|Ex^'=[wp(CT,x(-);CT) - Wp+i(CT,x(-);CT)]w(dx) < Ce "'^ \ 

Therefore, under the conditions of Corollary 7.1, and using the notation for the process, we 
have 

nClKlogK 

\{ho;Oixp~Xp+i))^\< / daE^g\wpia,^i-);a)-Wp+iia,^{-);a)\+OiK-^), (7.37) 

Jo 

where w is the solution to (7.17) with initial condition wo ~ V/io(x(ct)), assuming that the constant 
Ci in the upper limit of the integration is large enough. 

7.2.2 Restriction to the set Q and the decay estimates 

The complement of the set Qa.z includes the "bad" paths for which the level repulsion estimate 
in an average sense does not hold. However, the probability of QJ^ ^ is not very small, it is only 
a small negative power of K, see (7.25). This estimate would not be sufficient against the time 
integration of order CiKlogK in (7.37); we will have to use an — L°° decay property of (7.17) 
which we now derive. 

Denote the L^-norm of a vector u = {uj : j e /} by 

El"^n ■ (7-38) 

JG/ 

We have the following decay estimate for the equation (7.17): 
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Proposition 7.4 Consider the evolution equation 

dsu{s) = -A{s)u{s), u(s) e = M'^ (7.39) 
and fix a > 0. Suppose that for some constant b we have 

b 

and 



Bjk{s)>j—^. 0<s<a, j^fc, (7.40) 



W,{s)>^^, + 0<s<a. (7.41) 

Then for any 1 < p < q < oo we have the decay estimate 

||u(.s)||, < (.sfe)-(5-^'||u(0)||p, 0<s<a. (7.42) 

Proof. We consider only the case = 1, the general case follows from scaling. We follow the 
idea of Nash and start from the i^-idcntity 

ds\Hs)\\l = -2a{s)[n{s),n{s)]. (7.43) 

For each s we can extend u(s) : / M'^ to a function u(s) : on Z by defining Uj{s) = Mj(s) for 
IjI < K and Uj{s) = for j > \K\. Dropping the time argument, we have, by the estimates (7.40) 
and (7.41) with b=l, 

2a(u,u)> ^ (^!"^f >cHaH|H5H,-^ (7.44) 

with some positive constant, where, in the second step, we used the Gagliardo-Nirenberg inequality 
for the discrete operator \p\, see (B.l) in the Appendix. Thus we have 

a[u,u]>c||u||4||u||2-^ (7.45) 

and the energy inequality 

9.||u||2<-c||u||^!|u!|,-^<-c||ui|^||ui|^^ (7.46) 

using the Holder estimate ||u||2 < ||uj| j^''||u||4^^. Integrating this inequality from to s we get 

||u(s)|b<Cs-V2||u(0)||i, (7.47) 

and similarly we also have ||u(2.s)||2 < Cs~^''^||u(s)||i. Since the previous proof uses only the time 
independent lower bounds (7.40), (7.41), we can use duality in the time interval [s, 2s] to have 

\\n{2s)U<Cs-'/^n{s)h. 

Together with (7.47) we have 

||u(2s)||oo<Cs-i||u(0)||i. 
By interpolation, we have thus proved the Lemma. q 
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In the good set G (see (7.21)), the bounds (7.40) and (7.41) hold with b = cK^^' . Hence 
from the decay estimate (7.42), for any fixed a, Z, we can insert the other good set Qa.z into the 
expectation in (7.37). This is obvious since the contribution of its complement is bounded by 



CiKlogK 



daE"Q;^^2^;|it;p(cr,x(-);cr) - u;p+i(cr,x(-);cr) 



< CK^ 



CiKlogK 







da a-T^ W [g(^ |V,;/io(x(a))r+«) ^« Q%^, 



Act a i+e 



CxKXogK 



dcr a i+e 



{\ + d{xK{o))-^ ^-d[x_K[p))-^ 



-n3/2 



2/3 r 



(7.48) 



where we used the decay estimate (7.42) with q = c», p = 1 + ^ in the second line. In the third 
line we split the sum into two parts and used the bound 



|a,/io(x)|< |(Vy)'(x,)-(yy)'(x,)|< 



dixj)' 



(7.49) 



that follows from (6.42) (with ^ replaced by since y,y G TZL,K{£,^S/2,a/2)). Recall that d{x) is 
the distance to the boundary, see (6.44). For indices away from the boundary, \i\ < K ~ CK^ , we 
have |<i(a;i)| > v[mv{\i — \i + K\] on the set Q that guarantees the finitencss of the sum. For 
indices near the boundary we just used the worst term with i = ±K. We used a Holder inequality 
in the fourth line of (7.48) and computed the expectation by using stationarity and (6.63) in the 
last line. Hence we have proved the following proposition: 



Proposition 7.5 For any a > and for any path {x(-) : [0,(t] 
(7.17) with Hq given by (6.52). Suppose that 



p > 12^' + 6{Ci + C3 + 1)C 
holds with C4 defined in (7.25). Then for any fixed Z,p ^ I with p 7^ K , we have 



let w be the solution of 
(7.50) 



|(/io; 0{xp 

^p+i))ij\ — 



CiKloeK 



daE"Q,,za|^«p(a,x(.);a) -«;p+i(a,x(.);a)| +0(A-^/6)_ 

(7.51) 



7.3 Short time cutoff and finite speed of propagation 

The Holder continuity of the parabolic equation (7.17) emerges only after a certain time, thus for 
the small a regime in the integral (7.51) we need a different argument. Since we are interested in 
the Holder continuity around the middle of the interval / (note that \p\ < K^~^ in Theorem 6.5), 
and the initial condition V/iq is small in this region, a finite speed of propagation estimate will 
guarantee that Wp{a; x(0), cr) is small if a is not too large. 

From now on, we fix cr < CiKlogK, \Z\ < K/2 and a path x(-), and assume that x(-) e 
G n Qa,z- In particular, thanks to the definition of G and the regularity of the locations aj, the 
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time dependent coefficients Bij{s) and Wi{s) of the equation (7.17) satisfy (7.40) and (7.41) with 
b = K-^' . 

Let q(s) := w(s; x(-), a) which satisfies the equation (7.17) with initial data q(0) = Vft-o(x((j)), 
where /iq is defined in (6.52). We also recall the bound (7.49). We split the solution q = q™ +q°"* 
according to the support of the initial condition, i.e. q*" and q°"* are solutions to (7.17) with 
initial data 

g;"(0) := a>o(x(a))l(|j| < K^-'^^), g°"*(0) 9,/io(x(a))l(|j| > A'^-''-^) (7.52) 

with some small parameter ^5 > 0. From (7.49) and the rigidity bound given by t/, we have, on 
the set Q, that 

||q"(0)||i < A'^'-^N (7.53) 

From the — L°° decay estimate (7.42), we have 

/ da E'^g\q';\a) - q';'+,{<7)\ < K^'""^' / <j-^ da < K^^'-'^ (7.54) 

Jo Jo 

(Recall that the singularity a ^ can be cutoff as in (7.48)). Together with (7.51) and with the 
choice 

05 >P (7.55) 

and recalling p > 4^' from (7.50), we have 

\{ho;Oixp-Xp+i))^\< / daE''Q^,zg\q°/\'j)-q;i\ia)\+0{K-f/^). (7.56) 

Jo 

It remains to control (7°"* . In order to study this question, we consider the fundamental solution 
of the equation (7.17), i.e., for a fixed, the solution 

dsu''{s)^-A{s)u''{s), u'^{Q) = daj. (7.57) 

Note that the equation is the same as the one studied in (7.17), but the initial condition is different; 
in particular previously we studied this equation with a random initial data. By linear superposition 
we then have 

qTit)= dahoi^{a))u;it) (7.58) 

for any t > Q. 

We will assume that the coefficients of A satisfy for some fixed \Z\ < K/2 and pi > the 
bound 

sup sup / -J- E B,,ia~s)ds<CKP\ (7.59) 

and the similar bound if Z is replaced with K and —A'. Notice that on the set Q n Qa,z and with 
the choice 

Pi :=p + e', 

(7.59) follows from (7.22), since for the summands with \i — j\ > in (7.59) we can use Bij < 

The following lemma provides a finite speed of propagation estimate for the equation (7.57). 
This estimate is not optimal, but it is sufficient for our purpose. The proof will be given in the 
next section. 
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Lemma 7.6 [Finite Speed of Propagation Estimate] Fix a Cz I and a < CiKlogK. Denote 
by u(s) = u°(s) the solution to (7.57) with initial condition u°(0) = Saj. We assume that the 
coefficients of A satisfy 

W,(s)>^, Sy(s)>-^^, 0<s<a, (7.60) 

where di := inin{|z + ifj, |i — KW. Assume that (7.59) is satisfied for some fixed Z , \Z\ < K/2. 
Then we have the estimate for any s < a 

\u;{s)\ < ^^i^ (7.61) 

\p-a\ 

7.4 Proof of the Finite Speed of Propagation Estimate, Lemma 7.6 

Let 1 ^ £ ^ K he a. parameter to be specified later. Split the time dependent operator A = A{s) 
defined in (7.15) into a short range and a long range part, A = S + TZ, with 

(5v), :=- B,k{vk~vj) + W,v, (7.62) 

k : \j-k\<e 

and 

(7^v), :=- B,u{vu~v,). (7.63) 

k ■ \j-k\>t 

Note that S and TZ are time dependent. Denote by Us{si,S2) the semigroup associated with S 
from time Si to time S2, i.e. 

ds2Us{si.S2) ^ -S{s2)Us{si,S2) 

for any si < S2, and Us{si,si) = /; the notation U^{si, S2) is analogous. Then by the Duhamel 
formula 



u{s) ^Us{0,s)u„+ UA{s',s)n{s')Us{0,s')uods'. (7.64) 
Jo 

Notice that for £ ^ and for x(-) in the good set Q (see (7.21)), we have 

ll^-lli=E| E ., ^.|<g^-^l|vi|i, (7.65) 

\j\<K k:\j~k\>i ^ ' 

or more generally, 

\\nw\\p<cr^\M.p, i<p<oo, llTiviu <C£-2|iv||i. (7.66) 

Recall the decay estimate (7.42) for the semigroup C/4 that is applicable by (7.60). Hence we have, 
for s > 2, 

S 

||C/^(s',s)7^(s')t/5(0,s')uo|Lds' 

<A'«' / {s~s')-^\\n{s')Us{Q,s')n4^ds' <K^'£-\\ogs)\\n4i, 
Jo 
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where we used that Us is a contraction on L^. The non-integrable short time singularity for s' 
very close to s, \s — s'\ < K~'~^ , can be removed by using the — > L°° bound (7.42) with some 
p > 1, invoking a similar argument in (7.48). In this short time cutoff argument we used that Us 
is an L'P contraction for any 1 < p < 2 by interpolation, and that the rate of the — >■ L°° decay 
of Ua are given in (7.42). 



\vl{s) -Us{Q.3)xio\\oo < r\\ogs)K^' <Cr\\ogK)K^\ 



(7.67) 



where wc have used that x(-) is in the good set Q and that s < CiK log K . 

Wc now prove a cutoff estimate for the short range dynamics. Let r(s) := Us{0, s)uo and define 



/(.) = ^0,r|(.), cj^,=e\=-^\/' 



(7.68) 



with some parameter 6* > £ to be specified later. Recall that a is the location of the initial condition, 
Uq = 6a- In particular, /(O) = 1. 

Differentiating / and using Wj > 0, we have 

/'(s)-a.^0,r2(s)<2^0, r,is)Bk,{s){rk-r,){s) 

3 3 k:\j-k\<l 

= XI Bkj{s){rk - rj){s) [rj{s)(j)j - rk{s)(j)k\ 

\3-k\<i 

= XI Bkj{s){rk - rj){s)(j)j [rj - r^] (s) 

|j-fc|<£ 

+ X ^kj[,s){rk~rj)[s){^^- ^k\rk[,s). 
\j-k\<i 

In the second term we use Schwarz inequality and absorb the quadratic term in — rj into the 
first term that is negative. Assuming i < 9, we have <^^^ — 0fc]^ < for \j — k\ < £. Thus 



\j-k\<e 



k',j:\]-k'\<e 



From a Gromwall argument we have 



f{s) < cxp 



^ , , , ■ , I ■ lA^ I) 



k,r.\]~k\< 



/(O). 



(7.69) 



From the assumption (7.59) with M — K and any arbitrary Z . wc can bound the integration in 
the exponent by 

(7.70) 
(7.71) 



k,j:\j-k\<K 



Thus we have 



Xe'^"'''/^'(s) = /(s) <exp [e-^fK'+P^s + l)]fiO) < C, 
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provided that we choose 

e = eK'-P'+^'^/^V^Tl. (7.72) 

Now we choose 

^= b-a|X-«'-(''i+i)/2(s + l)-i/2 

so that e'P""!/^ > cxp {K^ ). Using this choice in (7.71) and (7.67) to estimate Up(s), we have thus 
proved that 

K{s)\ < £-\log K)Ki + Ce-^ < . p^^. (7.73) 

\p-a\ 

This concludes the proof of Lemma 7.6. 

7.5 Completing the proof of Theorem 6.5 

In this section we complete the proof of Theorem 6.5 assuming a discrete version of the De Giorgi- 
Nash-Moser Holder regularity estimate for parabolic equations with singular coefficients (Theo- 
rem 7.7 below). 

Notice that on the set Q n Qa-,z the conditions of Lemma 7.6 are satisfied, especially (7.59) 
follows from the definition (7.22) as we already remarked. Thus we can use (7.61) and (7.58) to 
estimate the short time integration regime in (7.56). Setting 

6*5 := mini — , — I, (7.74) 
^ I 2 ' lOOi' ^ ' 

we obtain, for any \Z\ < 2K^~^' and \p\ < K^~^' , 

da E"Q,.ze|gr*(a) - g°fi(a)| (7.75) 



<C I daWQ^^zG \daha{^{<j))\ul{a 

< Cif2«'+w+i/2+i+i-(i-e5) sup E^ Q^ zQ C ds V \daha{^{<j ^ s))\ 

< CK'^^'+P'+'''--^ sup E"^ Q^.zG f ds 

l<cr<Kl/4 Jo 



l<cr<Kl/4 



1 



d{xK{<J ~ s)) d{x-K{(J-s)) 



provided that 

^e + Pi + Cd<^- (7.76) 

In the third line above we used (7.61) together with |p — a| > ^K^^^^. This latter bound follows 
from \a\ > K^'^'^ and \p\ < K^~^' and from the choice O5 < £,* . In the fourth line we used (7.49) 
and that on the set G we have 



d{xj) id{xK) d{x^K) 



1 1 
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Moreover, in the last step we used (6.63). 

Wc now treat the integration regime a G [K^/-^,CiK\ogK] in (7.56) which requires a partial 
regularity result which we now formulate. We define the regular kernel by 



- (7.77) 



where we recall that p(0) is the equilibrium density at the center. Since x(-) G Q, we have 



sup sup 

s<CiK\ogK \^_j\yCKi' 



B,j{s)-B,j <K-^'b,j (7.78) 



for some sufficiently large C depending on £'(0). In particular, there is a constant C such that 

^ < Sy(.s) < (7.79) 



for any \i - j\ > CK^' and < .s < CiA'log A'. 

Theorem 7.7 (Parabolic partial regularity v^rith singular coefficients) Let u be a solution 
to (7.57) , where u = u" for any choice of a. Suppose that there exist ^o, Po small enough such that 
for some ^' < '"^'^ Pi l£ Pa the coefficients of A satisfy (7.60), (7.59), (7.79) and 

W,<—, if d,> K^^ . (7.80) 

di 

Let a G , Ci K log K\ he fixed, where C3 > is an arbitrary positive constant. Then for any 
< q' < f there exists q > depending only on ^q, Pq and q' so that for any \Z\ < K/2 

sup \uj{a)-Uj,{a)\<Ca~^~'^, (7.81) 

max(|j-Z|,|i'-Z|)<ffi-<!' 

where u = u° for any choice of a. The constant C in (7.81) depends only on the control parameters, 
C3,po,^Q, but it is uniform in a. 

Notice that this result is deterministic. We also remark that if we define the rescaled function 
v{j/K,t) tujit), then (7.81) can be interpreted as a type of Holder regularity of v on scale 
cri-?'/^— 1 < 1 at the point Z/K: 

\v{x, a) - v{y, a)\ < ct"' < |a; - yl''^"^ (7.82) 

for l/K < \x ~ y\ < cr^"'? /K and x,y near Z/K. The Holder exponent is thus at least c^q. 

We have stated the regularity for \Z\ < K/2 to avoid dealing with the boundary regularity. For 
the application in the proof of Theorem 6.5, we need the regularity only for Z = p (or any Z such 
that \Z — p\ is bounded). We will prove Theorem 7.7 in Sections 8. 

Assuming Theorem 7.7, we now complete the proof of Theorem 6.5. As we already remarked, 
the conditions of Theorem 7.7 are satisfied on the set Qa,z H Q with some pa,£,o small universal 
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constants. For any \p\ < ^ fixed, we choose Z = p (in fact, we could cfioose any Z witli 
\Z — p\ < C). Using (7.58) and (7.49), we liave, for the large time integration regime in (7.56), 

CiKlogK 

da Q^^zGl (^) - Ipti (^) I (7-83) 



JkU^ , ,.^_«, d(Xa(CT)) ' 



nCiKlogK _ <.! 

<CK^^' da a-^-'i E'^ zG ds V 



1 



d(xa(a — s)) 



< CK 



In the third line we used Theorem 7.7 with €3 = 1/4 and q' = 1/2, which guarantees that \Z — p\ < 
C < a^~'^ for any a > K^^^. The exponent q, corresponding to the choice q' = 1/2, is a small 
universal constant. In the last line of (7.83) we used a similar argument as in the last step of 

(7.75) . 

Finally, from (7.56), (7.75) and (7.83) and pi = p + ^' wc have 

\{ho;0{xp - Xp+i))^\ < CK^^'+P+'^'^-i" + 0{K'^) + 0{K-P/''). (7.84) 

For a given ^* > 0, recall that we defined ^5 := min{^*/2, j^} and we now choose 

. ( q 05\ • / 9 1 \ (7QK\ 
p:=mm<^ — , — > = mm — , — , >, 7.85 

llO 2 / llO' 4 ' 200/' ^ ' 

in particular (7.55) is satisfied. It is then clear that for any sufficiently small ^ all conditions in 

(7.76) and (7.50) on the exponents ^, ^' = (C2 + 1)C^ and pi = p + C can be simultaneously 
satisfied. Thus we can make the error term in (7.84) smaller than K'^^K~''/^ . With the choice of 
£ = p/6, where p is defined in (7.85), we proved Theorem 6.5. □ 

Although the choices of parameters seem to be complicated, the underlying mechanism is that 
there is a positive exponent q in (7.81), which to a great degree is an universal constant. This 
exponent provides an extra smallness factor in addition to the natural size of Uj{a), which is a~^ 
from the L°° decay. As (7.82) indicates, this gain comes from a Holder regularity on the 

relevant scale. The parameters ^, ^' and ^* can be chosen arbitrarily small (without affecting the 
value of q). These parameters govern the cutoff levels in the regularization of the coefficients of 
A. There are other minor considerations due to an additional cutoff for small time where we 
have to use a finite speed estimate. But the arguments for this part are of technical nature and 
most estimates are not optimized. We just worked out estimates sufficient for the purpose of 
proving Theorem 6.5. The choices of exponents related to the various cutoffs do not have intrinsic 
meanings. 

As a guide to the reader, our choice of parameters, roughly speaking, arc given by the following 
rule: We first fix a small parameter Then we choose the cutoff parameter 9^ to be slightly 
smaller than (7.74). The exponent p in (7.22) has a lower bound by ^ and ^' in (7.50). On the 
other hand, p will affect the cutoff bound and so we have the condition p < 65 (i.e., (7.55)). So 
we choose p ^ and make ^, ^' very small so that the lower bound requirement on p is satisfied. 
Finally, if the parameter ^* < q/10, we can use the gain from the Holder continuity to compensate 
all the errors which depend only on 
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8 A discrete De Georgi-Nash-Moser estimate 



In this section we prove Theorem 7.7, which is a Holder continuity estimate for the parabohc 
evolution equation (7.57). Our equation is of the type considered in [10], but it is discrete and 
in a finite interval. The key difference, however, is that the coefficient By- = (xj — Xj)~'^ in the 
elliptic part of (7.57) can be singular, while [10] assumed boundedness. In fact, by extending the 
reasoning of Ben Arous and Bourgade [3], the minimal gap min,i(a;i+i — Xi) for GOE is typically 
of order N~^^^ in the microscopic coordinates we are using now. The only control we have for the 
singular behavior of Bij is the estimate (7.59). This estimate essentially says that the space-time 
maximum function of BiA+i (t) at a fixed space-time point (Z, 0) is bounded by K^. Our main task 
is to show that this condition is sufficient for proving Holder continuity at the same point. Our 
strategy follows closely the approach of Caffarelli-Chan-Vasseur [10] and the key cutoff functions 
(8.1, 8.30) are also the same as in [10]. The main new feature of our argument is the derivation of 
the local energy estimate. Theorem 8.1, for parabolic equation with singular coefficients satisfying 
(7.59) and (7.60). The proof of Theorem 8.1 will proceed in two steps. We first use the condition 
(7.59) and the argument of the energy estimate in [10] to provide a bound in L^(L^(Z)) on the 
solution to (7.57) (Lemma 8.2). Using this estimate we run the argument again to improve it to an 
L°° estimate in space to obtain Theorem 8.1. Besides this proof, the derivation of the second De 
Giorgi estimate (Lemma 8.5) is also adjusted to the weaker condition (7.59). After that, starting 
from Section 8.2.2, the argument of [10] takes over and the proof presented here is the same as in 
[10]. Some more explanation will be given along the proofs. 

We warn the reader that the notations of various constants in this section will follow [10] as much 
as possible for the sake of easy comparison with the paper [10]. The conventions of these constants 
will differ from the ones in the previous sections, in particular, we will restate all conditions. 

8.1 Local dissipation estimate (first De Giorgi lemma) 

Fix a large integer number M and a center Z E I with M < Z < K — M. Let B = B^ z = {* G 
I \i — Z\ < M} be the interval of length 2M about Z, and let it be equipped with the counting 
measure. 

For any t> Q define the function 



(this corresponds to the case s = 1, = 1 with the notation of [10]). Note that this also has size 
of order i and localizes on a scale M . Let 



Here t will thus play the role of the typical size of u — ijj. One could scale out H completely, 
but we keep it in. For any real number a, we will use the notation a+ = max(a, 0) > and 
a_ — min(a, 0) < 0, in particular a = a_|_ + a_. In this section, we will use the convention 
that the time variable is in some interval [— Cif,0] where Ck is some constant depending on K. 
This convention is widely used for parabolic equations and in particular in [10]. Later on in our 
application, we will need to make an obvious shift in time. 

Theorem 8.1 (Parabolic energy estimate v^rith singular coefficients) There exists a small 
positive constant £o with the following properties. Let u he the solution to (7.57) and let M be 




(8.1) 
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defined by M"^ = K with some •& >\. Set a < CiK log K . We assume the following conditions on 
the matrix elements of A = B + W : 

i) For a fixed \Z\ < K/2 and p > assume that 

sup sup r E E B,,{a~s)ds<CK'>. (8.2) 

a) For some £^ > Q we have 



~dr 

where di :— min{|i + ifj, |i — /v |} 



W,{s)>——, B,^.(s)>— — , 0<s<a, (8.3) 



Hi) 

1 „ , . C 



<B,j{s)<——, (8.4) 



for any \i — j\ > CK^ and < s < a. 
Suppose 

sup ||u(t)||oo <^M^^ (8.5) 
te[-2J\/,o] 

for some Ci sufficiently small such that 

$ := 3ci + (3e + p)i? < 1/10. (8.6) 

Then for any K > Kq with Kq depending only on 'd and $ the following statement holds. For any 
Z and i as above if 

i /°^^dt5](^.(0-^.)+]'^'<eo^, ^^ = 4^-^''\ (8.7) 
(sq was given in the beginning of the this theorem) then for any \i\ < K we have 

e 

sup u,(<)<-+V, (8.8) 

t£[-ALQ] ^ 

We remark that the constant Sq is independent of any other parameter, in particular it is indepen- 
dent of £, Z, M,ci,^, pj-d. 

The bound (8.8) is essentially an — ?> L°° estimate, which is the basic input for applying the 
method in [10]. Notice that (8.2), (8.3) and (8.4) are the same as (7.59) (7.60) and (7.79) but gi 
and ^' are replaced by g and ^ for simplicity of notations. 

Proof. We first prove a weaker version of this result, namely, the following lemma: 

Lemma 8.2 Under the assumptions and notations of Theorem 8.1 we have 

sup y^{u,{t)-i^l'%<CM^l\ (8.9) 

tG[-A/,0] j 

where $ is defined in (8.6). 
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Proof. Assume for notatioiial simplicity that Z = (this requires relabelhng the configuration 
space). By direct computation, we have 

^t\YM^ - iiW = - Y.^u, - i^D+B^A^^ - ^3) - E("' - ^')+^'"'' (8-10) 

i ij i 

Recall that Bij depends on time, but we will omit this from the notation. Since Wi > 0, the last 
term can be bounded by 

i i 

In the first term on the right hand side of (8.10) we can symmetrize and then add and subtract 
■tp^ to u we get 

-^{uz - ^i)+Bij{ui - Uj) = - b[{u - tl/)+,u] 

ij 

= - b[{u - /)+, {u ~ v')+] - biiu ~ v')+, - V-')-] - Hiu - v/)+, V' 

Since Bij > and [a+ — 6_|_][a_ — > for any real numbers a,b, for the cross-term we have 
b[{u — ip^)+, {u — i^^)-] > 0. Thus the last equation is bounded by 

< -b[{u - {u - V')+] - b[{u - ^']. (8.11) 

Using the definition of o (7.16), we have thus proved that 

EK - < -a[{u V')+, (u - b[{u - V-']- 

i 

Decompose the error term into 

b[(u-i^^)+,V^] = r!i + r!2 + 03, 

^^1^=^ E Sy■[^'f-V',']((^^^-^^)+-(%-^',')+) 

\i-j\>M 

and ft2 and ft^ are defined in the same way except that the summation is restricted to CK^ < 
|i - j\ < M for and \i - j\ < CK^ for fig, where C is the constant from (7.78) and (7.79). We 
have 



\i.i-i.]\< ^r^:' (8.12) 



The last term is bounded by £AI — Recall from (7.79) that in the regime \i — j\ > M 

we have Bij < jjrjp- Together with M > K^, we have 

\n,\<m-'/' E ^Twt("»-^')+ + ("^-^.')+]^]^EK-^')+ 

\t-j\>M I i 
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For il2 1 by symmetry of Bij we can rewrite it as 

dKi<\i-j\<M,lf>f<i:'- 



4 



Cif«<|i-j|<A/ 
CKi<\i-j\<M 

The first term is ^b[(it — ^ "0^)+] ^-nd can be absorbed in the first term on the r.h.s. of 

(8.f f ). By (8.f2) and (7.79), the second term is bounded by 

4 E S.,|V'f-0jl'x(".-V^f>O)<C£2Af-i5]x(^.-V'f >0), (8.f3) 

CK<^<\i-j\<M * 

i.e. 

using that b < a. 

A similar estimate is performed for fis, but in the corresponding last term we use that 

\i>\-i^'^<C{llMf 

for |« — j| < ■ Thus we have 

|i-jl<Cife 

< Cif 2« (f/M)2 ^ - > 0)B,, 

\i-j\<CKe 

^ - V4 > 0)[i3M+i (8.14) 

i 

Here we just overestimated sums by CK^. The conclusion of the energy estimate is 

i 

CP . __ 
+ M E("» - ^')+ + E ^("^ " '^'^ > 0) + "4 (8.15) 
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(8.16) 



Due to the assumption (8.5), we can assume that the summation of i is restricted to \i\ < M^'^^'^'^. 
By (8.2), we have 



M1-3C1 



2 /-t 



ds 



^ [i3,.,+i+i3,,,_i](s)<(t-T)i^3€+P£2^/3ci-i. (8.17) 

|i|<A/i + 3':i 

For a fixed A/, let Tk = -1/(1 + 2-*^), 4 = 1(1-2"'^') |, = [rfc,0] x {1,2,... A'}. Define 



1 1 



[{u-i;'^)+,{u-i,''')+]{s)ds 



The prefactors 1/(M^|) just make the size of U order one. 

Integrating (8.15) from r to t with Tk-i <t <Tk <t <0,we have from (8.17) 

- ^Ml W + /* a[{u V^)+, (« - V>)+](s)ds 



+ C{t - t)A'3«+''^2^/3=i-1 



^ 5](u, - ^^)+(s) + ^ 5] x(". - i^^ > 0){s) 



M 



ds 



(8.18) 



Taking the average over r e [Tfc-i, Ta;_i + 2-'=-iA/r] = [Tfc_i, - 2"''"^A/] and using that in this 
regime t ~ t ^ 2~^~^M, we have 

- ijf^]l{t) + f a[{u - ^A)+, (^, _ V'^'=)+](s)ds 
2fc+i rTk-2-''-'^M 



< c- 



M 



Tk- 



Y,[u^-i^^]l{s)ds 



C 



■2 
7\/ 



|i ^(«, - ^,f'=)+(s) + ^ - > 



ds + CK^^+Pi^M^"^ 



Dividing through by and taking supremum over t e [Tfe, 0], for fc > 1 we have 



ofc+l /-o r 1 1 



(Ui - 7/'^'=)+ + X{U^ - tpi" > 0) (s)ds 



M 



(8.19) 



The first three integrands have the same scahng dimensions as jP'. One key idea is to 
estimate these terms in terms of the L'^ norm of u and then using the Sobolev inequahty. It is 
elementary to check these three integrands can be bounded by the L** norm of (u — -i/'^'°)+i by using 
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that if u, > V'-N then m - i/^f'-' > 4 - 4-i = > 2^^''+^k: 

- <E(^'' - V'-'")+x(«. - > 2-^'+'k) (8.20) 

i i 

i 

- # > 0) <(2^-+^)VE(-^ - ^'r')' 

i i 
i i 

Recall the Sobolev inequality (7.45) ( //) < ( Ei /i^) /]^^^- Thus we have proved 

that for fc > 2 

Uk <(2'=+^)^^ d.[EK - a[(« - {u - ^'^-^)+]{s) 

<C^Ul_^ + fc > 2. (8.21) 

recalling that K ^ M'^ and $ 3ci + (3^ + p)z?. 

For A; = 1, we estimate the integrands in (8.19) by £^-norms. We have the following general 
estimates for any i' < I" 



Y^(u, - V.r)+ < E("^: - 4)^X(U^ - > ^" - < ^77^77 E^"' " )+ (^.22) 

i % i 

E - € > 0) <7^;7:^ E("« - 

i ^ ^ i 

We use (8.22) with I" = ti and ^' = in (8.19), this implies that 



Without loss of generality, we assume that Ci > 2, where Ci is the constant in (8.21). Now choose 
the universal constant eo in (8.7) so small and M big enough so that this last inequality implies 

1 

Ui < 



loocf ■ 

One can check the recurrence relation (8.21) implies that 

Uk < (8.23) 

for any fc such that C^'^ < M~^+*. Choosing fc sufficiently large, noticing that Uk is decreasing 
in fc (as Tk and £k increase), we find from the i!^-norm term in Uk that (8.9) holds. This proves 
Lemma 8.2. q 
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We now repeat the proof of Lemma 8.2 but with £k, fc > 1, replaced by 



2e, 



4 = -(1-2-^-^) 
5 



(8.24) 



and with integration on a smaUer time intervaL Choose Mi such that Af* ^ Mi <C M^ *, and 
set fk := -Mi{l + 2-^). Define 



a[(ti-7A^'=)+,(w-^A^^)+](s)d 



The proof is unchanged up to (8.17), the integral of 



M2 



Y^xMs) - > 0)[i?,,,;+i(s) + B.,,-i(s)] 



1+* 



is still estimated by 

j tli{s)ds < C{t - T)K^^+PpM-^+^''^ < Cit - T)e^M-^ 
Similarly to (8.18), we integrate (8.15) (with ^replacing g) from t to t with fk-i <T<fk<t<0 
J2[u, - V|li(i) + /* a[{u - /^)+, {u /'=)+] (s)ds 

^ J2{u, - 4>')+{s) + ^ 5] x(". - > 0)(.) 



<5^[u,-V'M^(r) + C 

i 

+ c{t^T)eM-^+^. 



ds 
(8.25) 



Taking the average over r G [Tk-i,Tk-i + 2~''~M/i] = [Tfe_i,Tfe - 2'''-^Mi] and using that i 
this regime t — t ^ 2~^~^Mi, we have 

- ^f'']^(t) + r a[(ii - /'^)+, {u - VA)+](s)ds 

2fc+l /•Tfc-2-'=-iM: 



< c- 



Ml 



c 



A/ 



i 



ds + CfMiM-^+'^ 



Dividing through by A/^ and taking supremum over t G [7^,0], for fc > 1 and using Afi < Af, we 
have 



2fe+i /-o 2^ 



Uk <c 



K - + y{u, ~ 4")+ + x(«. - 4" > 0) (s)ds 



CMiM' 



-2+* 



(8.26) 
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Using the bounds (8.20) and the Sobolev inequahty, we get, instead of (8.21) for any k > 2 



Uk < 



M ^ 

<Cf [/? 1 + Cil/iAf^2+*, k>2. 

- 1 Ml '--1 ^ ' - 



(8.27) 



This new recurrence inequahty has the solution 

Uko < CA/iA/"2+* (8.28) 
for a sufficiently large /cqi as long as the recursion can be started, i.e. if we knew 

Ml 



Ui < 



M 



29) 



For k — 1 the estimate (8.26) together with (8.22) (with £i replacing ii) becomes 



Ui < 



C 



<- 



MM^ 
C 



1 J-2Mi 



- + - 4')+ + X{u^ - ^-^ > 0) (s)ds + CMiM 

1 ^1 



<- 



-MM- 

.CM" 
M 



-2+a 



-2+* 



-2A/i 



In the second step we used (8.22) with £" = ii and £' = £/3 noting that £i = ^£ > \£. In the last 
step we used (8.9) and 2Mi < M. Thus (8.29) is satisfied if Mi > M*. 
Finally, (8.28) implies 

sup y^{u,{t) - < C£^MiM-^+^ = o{£^) 

te[-Mi,o] j 



since Ah < A/^"*, i.e. 



sup snp{ui(t) - ■0^^^)+ < o{£) 

t£[~MuO] i 



The estimate (8.8) now follows if we repeat the same argument for each time interval of length Mi 
within [— Af, 0]. This completes the proof of Theorem 8.1. q 



For any Z, M and £ define 
We have the following corollary 



(ALZ.e) 



i~ z 


1/4 X 






M 





.30) 
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Corollary 8.3 There exists a constant 5 > Q (depending only on from Theorem 8.1) such that 
for any K > Koi'^, ^) if for some M, Z and £ a solution u satisfies (8.5), 



sup max Ui {t) < i + ipi 

tG[-2A/,0] ' 



and 



the 



[ #\ie[Z- 2M, Z + 2M] ; u,{t) > o]dt < S, 



sup max Ujit) < -. 



(8.31) 



(8.32) 



.33) 



Proof. Fix (^ojio) £ [^-^-^,0] x |— M, A/], our goal is to show that ■tiio(to) < ^/2. Consider the 
interval /q = Jxq — M/R, xq + M/i?] of length 2M/R about xq, where i? is a large constant. We 
will apply Theorem 8.1 to the interval Iq, i.e. we define 



i - io 1/2 



M' 



1 



M' := M/R. 



Since V'io = ^^ the estimate ^^^(to) < £/2 would follow from the conclusion of Theorem 8.1. We 
just have to check that 



1 



{M'Y 



2M' 



1/2 



.34) 



if 5 is sufficiently small (depending only on eo). Notice that if i? is a sufficiently large universal 
constant, then 

^ + V'i < V'i, for \i\ > 2M. (8.35) 
This follows from the simple inequality 



1 



i 1/4 



M 



-1 < 



I - iQ 



M/R 



1/2 



- 1 



whenever \io\ < M and |i| > 2M. Thus we have 



/" y^{u,it) ~ VO+di = /" y\ Mt) - i^^)ldt < M^S max {£ + i,,f < C2M^eH. 

with some universal C2. In the first equality we used that for \i\ > 2M the condition (8.31) and 
(8.35) imply that Ui{t) < ipi. In the second inequality we used that {ui — ipi)+ < (wi)+ < i + 'ipi 
and (8.32). The last inequality is a simple arithmetic estimate. Thus (8.34) follows if 5 is so small 
that 

C2i?2<5 < el 



This completes the proof of the Corollary. 

Remark 8.4 Suppose that in addition to (8.3), (8.2) and (8.4), the bound 

< — , tf d, > 



□ 



(8.36) 
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analogous to (7.80) is also satisfied. For any constant \A\ < lM^+''\ let V — u — A. Then v 
satisfies 

dsVr{s) = -{A{s)v{s)), - AW^ (8.37) 

We can apply the proof for Theorem 8.1 to v and view WiA as an error term. Then instead of 
(8.8), we will have 



sup f. 
te[-M,o] 



(8.38) 



provided that 7? > 3. The last error term IK ^1"^ is much smaller than the Holder continuity bound 
we aim for Theorem 7.7 and can be neglected. Similar remark applies to Corollary 8.3. 



8.2 Local decrease of oscillation: Second De Giorgi lemma 

As before, we are given three parameters, M, Z, £. The eenter Z ean be set zero for notational sim- 
phcity (this requires relabelhng the configuration spaee). Define a nonpositive function supported 
on \i\ < 3M by 



F, = F 



max — 1, min I 0, 





i- z 




(0, 






M 



Notice that < F < and i^; = if |i - Z| > 3A/. For any A > smaU, set 
i,, = ^'^.^.^) ^ [ ( 1^ -\-^y' -l for \i-Z\> 



(8.39) 



and V'i = otherwise. We also define three cutoffs, all depending on all four parameters, M, Z, ^, A: 



(0) 



(1) 



e + ij., + XF, 



Notice that 



(8.40) 

and when \i\ > 3M all inequalities becomes equalities. Notice that ipf^"^ = if |i — < 2y/2M for 
k = 0,1,2. 

Lemma 8.5 Let = K and assume that (8.2), (8.3) and (8.4) hold with some p and ^ so that 
(8.6) holds. Recall the constant 5 from Corollary 8.3 which depends only on Sq from Theorem 8.1. 
Depending only on this 6, there exists /i > 0, 7 > and A G (0, 1/8), such that whenever a solution 
u(t) to (7.57) satisfies 

u,{t)<i + ^,, <e[-3A'/,0], Vi (8.41) 
I 3M ^[Z-M,Z + M] : u,{t) < }di > /i, (8.42) 



and 



.43) 
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the 



-3M 



(8.44) 



If in addition (8.36) is satisfied and v is a solution (8.37) satisfying (8.41), (8.42) and (8.43), then 
(8.44) holds. 

This Lemma asserts that whenever the function u increases from tp^^^ to ip'^'^^ in time of order 
M, then there is a time interval of order M so that u is between tp^^'> and ip'-^\ The rest of this 
section is devoted to the proof of this Lemma. 

Proof. We can fix /i < 1/8. The parameter S may have to be reduced a bit along the proof. 
We are looking for a sufficiently small A so that there will be a positive 7. The key ingredient of 
the proof is an energy inequality (8.48) including a new dissipation term which was dropped in the 
previous section. Wc first derive this inequality. 



8.2.1 Dissipation with the good term 

Let -3A/ < Ti < < 0. We use (8.11) (with cutoff (^(i)) but we keep the good term 



T2 



b[(u(t)-(^(i))+,(u(0-¥'''')_]dt- / ' b[(M(<)-(^(i))+0,(^(i)]dt 



Ti 



(8.45) 



Notice that we inserted a characteristic function Oi = l(|i| < 3Af) using the fact that (8.41) and 
(8.40) imply Ui < v'^^'' for |i| > 3M, i.e. Ui — lyS^^-* = {ui — Lp^^'^)9. The last error term will be 
Schwarzed and partly absorbed in the quadratic term in the left. Here we use 



b{fe,g) = Y.^fA - fj0j)B,,{g, - g,) = ^(/,0, - fje,){e, + e, - 0,e,)B,j{g, 



9], 



so 



Hfe,g)\ < - M'B,, + ^ e,B.,,{g, - g,f 



Ti 



In the other term we use again a Schwarz to separate the two summands in ip^^\ We are left with 
to control 

>^''J2(F,~FJfB,,+J2ii^^-i^JfB,Je, it)dt (8.46) 

Since \Fi — Fj\ < C£M~^\i — j\ and supported on \ j\ < 3M, by splitting the summation to 
K ~ il ^ regime and its complement, we can bound the first term by 

r \^Y.^F,-F,fB,,{t)dt<\HHr^ r N-jf^».W 



T2 



<X^fM~^K^^ J2 B,^,+iit) + CX^e^M-^ 



j|,|i|<3A/ 

T2 



1 |i|<3M 



E 



1 |i|,|3|<3J\/ 



I^-Jl 
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where we have used Bij < in the first regime and (8.4) in the other regime. By (8.2), we 

can bound the last Une by 

X^fK^^+P + \^i^{T2 - Ti) < CX^fiT2 - Ti). 

For the second term in (8.46), we use that ipiOi = and the supports of 9i and tpj are separated 
by a distance of order M ^ K^. Thus we can use (8.4) to estimate tire kernel: 

•^^1 i,j "^^1 |i|<3A/ |j|>A/A-4 ' 

<M{T2-Ti) <^'A2(r2-Ti), 

\j\>M\'-^ 

where we have used tpj ~ ^(j'/M)-^/^ for large j. 

Inserting these two error estimates into (8.45), wc have 

\ Y}um - b[{u{t) - v:,(i))+, {u{t) - ^(i))+]dt (8.47) 



t=Ti 



Define 



< - / b[{u{t) ~ ¥^^1))+, {u{t) - ^(i))_]di + Ce\\T2 - Ti). 



Since \ui{t) — ^p^p \ < \9i, we have, for all t, 

H{t) < CX^M. 

Hence we have 

H{T2) + / ' b[(u(i) - (p(^))+, {u{t) - v?(i))_]dt < H{Ti) + Cf\^{T2 - Ti) < Ct^X^M (8.48) 

for any -3Af < Ti < T2 < 0. Notice that b(/+, /_) > for any function /. 

We have completed our main task in using the estimate (8.2) to control local singularity of Bij. 
The rest of the proof for Theorem 7.7 is essentially the same as in [10]. For the convenience of the 
reader, we present the detailed proof. Besides providing some details and tracking the dependence 
of various constants and exponents for our application, wc follow the equations in [10] very closely. 

8.2.2 Time slices when the good term helps 

Let E be the time that u{T) is substantially below (p'^'^\ i.e.. 



S := 



jre (-3A/,-2Af) : #{[*[ < M : u,(T) < ^|"^} > \^iM'^. 
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> -cM-' [ Y^Mt) - W - ^fVdt (8.50) 



We have from (8.42) that 

|S| > ^M/x (8.49) 

By (8.48), we have 

CX'fM> f fa[(u(i)-(^(i))+,(M(t)-(^(i))_]dt 



where we have used that Ui{t) — ip\^^ is supported on \i\ < 3Af and 

B,,{t)>cM-\ |z|,|j|<3M, (8.51) 

with some positive constant c (this foUows from (8.4) and M ^ K^). For t G E the number of j's 
in |j| < M such that Uj{t) < ip'^^ is at least jfJ-M; for those j's we have 

-(«,(0-^f')->^«-^f >1-A2>i. 
Thus we can bound (8.50) by 

where we have used that {ui(t) — ip[^^)+ < A. 
Altogether we have proved 

if A is sufficiently small (recall that /.i is fixed) . Thus there exists a subset 8 C S such that 

|e| < A^/^M 

and we have 

i 

Choosing A small and recalling (8.49) we see that 

^(«,(t)-^f))2^ <A3-i£2^ (8.52) 

i 

holds on a set of times t's in [— 3A/, — 2A/] of measure at least Mfj,/8. 
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8.2.3 Finding the intermediate set 

Suppose that (8.43) is not satisfied, then there is a Tq G (^2A/, 0) such that 

#{« : (w,(ro) - ^)+ > O} > ]^M5 (8.53) 
and choose a Ti G E (tlien Ti < Tq) such that 

(such Ti exists by the conclusion of the previous section, (8.52)). 
We also have 

H{T,) = J2iu,in) > J2i^P - ^['^)ll{iu,iTo) ))+ > 0) 

i i 

> Y,f{\ - \^fF^l{{u.,[T,) - ))+ > 0) > Cf^^HHI (8.55) 

i 

with some positive constant Cp. This follows from (8.53); notice first that the set (8.53) must lie 
in [-3M,3M] (see (8.40) and (8.41)), and even if the whole set (8.53) is near the "corner" (i.e. 
close to i ~ ±3M), still the sum of these FiS is of order S^M since Fi is linear near the endpoints 
i = ±3A'f . 

Choose now A small enough (depending on the fixed 5) s.t. 

- 16 

Since H{T) is continuous and it goes from a small value H{Ti) < jgCpX^i'^S^ M to a large value 
H{To) > ^CpX'^i'^S^M, the set of intermediate times 

D:=\te (Ti,To) : ^CpX^i^S^M < Hit) < ^CfX^I^S^m] 
L 16 4 J 

is non-empty. 

Lemma 8.6 D contains an interval of size at least cS^M with some positive constant c > 0. 
Moreover, for any t € D, we have 

#{i : iff^ <u,it)} < ^SM. (8.56) 

Proof. By continuity, there is an intermediate time T' e -D C [Ti,To] such that H{T') = 
|CfA^£^(5'^M. We can assume that T' is the largest such time, i.e. 

H{t) > IcpX^fS^M for any t e [T',Tq] D D. (8.57) 
8 

Let T" = T' + cS^M with a small c> 0. We claim that [T\ T"] C D. For any t e [T' , T"] we 
can use (8.48): 

H{t) < H{T') + CfX^it - T') < -CpX^fS^M + Ccf-X^5^M < -CpX^i^S^M (8.58) 
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if c is sufficiently small. This means that as t runs through [T',T"], H{t) has not reached 
\Cf>?(-'^5^M, in particular [T ,T"] C (Ti,To) since H{To) is already above this threshold. Com- 
bining then (8.58) with (8.57), wc get [T',T"] C D. This proves the first statement of the lemma. 

For the second statement, wc argue by contradiction. Suppose we have : ipf'^ < Wi(T)} > 
i(5M for some t E D. Going through the estimate (8.55) for Tq = r, we would get H{t) > 

Cf^£^6^M, but this contradicts to r e D. This completes the proof of the lemma. q 
Define the exceptional set J" C -D of times where u is below t^^^-*, i.e. 

T:={teD : < 2M : u,{t) - ipf < 0} > ^iM^ 

This set is very small, since from (8.48) we have 

J-iM 



> 



J2 {um-^Yi+B,m^At)-vf^)-<^t 

\i\,\j\<-iM 



> -cM' 



\i\,\j\<3M 

^^^^ / E Kw-^r^di 



i|<3A/ 

where we restricted the time integration to T in the first step, then we used (8.51) in the second 
step. In the third step we used that whenever Uj{t) — Lp^p < (see the definition of J-"), then 

-(u,(i)-^f )_>^(1-A)>|. 

By (8.41), {ui{t) — < and {ui{t) — = if \i\ > 3A/. Hence we can continue the 

above estimate 

'^'>''''" ^ 2I7A X ?(".(') - ^S")i<i' - ihx L ^ f^^'^*'"!-"! 

Here we used that T C D and in D we have a lower bound on H{t). The conclusion is that 

with some fixed constant C. Using that \D\ > c6^M from Lemma 8.6, we thus have 

if A is sufficiently small (like A < c6^^). 

This means that |Z) \ > ^6^M. Now we claim that for ^ e D \ we have 

M 



A{t) := #{z : < < ^f^} > -. (8.59) 



69 



This is because t ^ T guarantees that the lower bound ^)\ < Ui{i) is violated not more than 
/iM < M/4 times among the indiees |i| < 2M. By (8.56), the upper bound Ui{t) < "y^i is violated 
not more than \5M < M/4 times. 
Finally, integrating (8.59) gives 

/° #f* : < u,{t) < ipT^}dt = [° Ait)dt > ^\D \ J-| > cSHd^. 

J-3M ^ ' J-3M ^ 

with some small c > 0, whieh implies (8.44). This proves Lemma 8.5. q 



8.3 Oscillation estimate 

For e > 0, A > define 



= i>, 



= 1 



i-Z 



M 



- 1 



for 



\i-Z\> MX 



+ 



and ^/^i = otherwise. This is essentially the same function as -i/!, but the growth at infinity is not 
but only lil"^. Define for any set S and any function / the oscillation Oscsf := sup^ / — infg /. 

Lemma 8.7 Let Z e p/, /v - Mj and £ > 0. Let = K and assume that (8.2), (8.3), (8.36) 
and (8.4) hold with some p and ^ so that (8.6) holds. Let X > be given in Lemma 8.5 and use this 
A in the definition ofip. Then there exist e > and A* > such that for any solution satisfying 



we have 



£ - < u,{t) < e + i'^, Vi, t e [-3Af,0], with ^ i^'^^^'^''^'^'^^ (8.60) 



OscQ,,u < (2 - X*)£, Qm [-A/,0] x [-M,M]. (8.61) 



The constants e and A* depends only on fi, 7, A from Lemma 8.5, 6 from Corollary 8.3 and Eq from 
Lemma 8.1. All these constants eventually depend only on Eq provided that K > A'o(i^, with 3> 
satisfies (8.6). 

Proof. This proof closely follows the original proof of [10]. We can set Z = 0. Without loss of 
generality, we can assume 



1 

M2 



-2M 
'3M 



#{* : \i\<M, u,it) < ipf''>}dt > ^^ 



(Otherwise we can take —u). Let x = [\i/M\ — A '^Y ^-i^d a = l/(4£) ^ 1. Then for a; > 1, 

= i(x - 1) < £^ '- = Aeij, 

a 

Choose fco such that A^*"'" = e. Then we have 

A-^'^V, < V-i 

For e sufficiently small so that fco = 1007^^ where 7 is from Lemma 8.5. 



(8.62) 



(8.63) 



..64) 
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For k < ko, define by the following equation: 



/fc+i) 1 ^y(fc) u(o)=M. (8.65) 



By (8.64), the upper bound in (8.60) becomes 



^(0 <i+ ^V'i < ^ + V'* t& [-3M, 0], Vi, k<ka. (8.66) 



(fc) (k) 

Notice that for any \i\ < M the sequence u) in k is decreasing, in particular u\ < £. This follows 



from induction and from (8.60). Moreover, from (8.62) we have 

-2M 
-3M 



1 p-IM 



since the set in question is increasing in k. 

We wish to apply Lemma 8.5 to u^'^^ which is an affine transformation of a solution u. As 
we note in Remark 8.4, subtracting a constant from a solution no longer yields a solution to the 
equation (7.57). However, since Wi is very smaU away from the edges by assumption (8.36), all 
estimates carry through with a negligible errors. Hence we can apply Lemma 8.5 to u^'^^ to get 
that if 

^ J[^ *{' ■■ "^^''W > ^r'}d^ > S (8.68) 
^ /°^^ #{^ : f^'^ < n^\t) < )}dt > 7. (8.69) 



then 
Therefore 

^ J[^*{^ ■■ nf\t) > ^f)}di < ^ /°^^#{^ ■■ uf\t) > ^f^yt-^ (8.70) 

Notice that by (8.66) and Ei = Qii \i\ > 3M, the inequality ^(t) > (pf^ can hold only if \i\ < 3M 
for any k < kg. Assuming |i| < 3Af and u['^\t) > (pf\ we have 

^iut'^-e) + e = uf'\t)>^[°^ (8.71) 

Since \i\ < 3M < X'^^M, we have, together with (8.40), that 

uf-^^ >\^{i^, + F,)+i>^f\ (8.72) 
Therefore, we can bound the last integral in (8.70) by 



1 



M2 



#h : uf(t)>^f'}di 
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We have thus proved that 

< !,,,*{' ■ l"l - ™' > vf'}<it-r (8-74) 

Iterating this estimate k times, we would get 

-L r : > )}dt < ^ /° #{* : \^\ < 3M, v<f^(t) > 'jdt - kj 

J — 3-A/ ^ — 3-A/ 

which becomes negative if, say, A;7 > 20. 

Thus there is a fc < fco such that (8.68) is violated, i.e., 

^ /"^^ : ^(i) > ^r}di < <5. (8.75) 

Our goal is to apply Corollary 8.3 to We first check the assumptions of this Corollary. 

From (8.66) and the definitions of ■0 and in (8.30, 8.39), we have 

^(fe+i)^^-) < £ + < ^ + i € [-3M, 0], Vi, fc + 1 < /so- (8.76) 

This verifies the first assumption (8.31) of Corollary 8.3. Since ^pf^ = for |i| < 2M, we have 

: NI<2M, .,r)(0>0}di 

where we have used (8.73) in the last inequality. This verifies the second assumption (8.32) of 
Corollary 8.3. Thus from Corollary 8.3 we get 

u[''+^\t) < e/2, t e [-A/,0], \i\ < M. 

Together with (8.65), we have proved 

u,{t) - £ = A2(^+i)(iif +i)(t) ~i)< -A2(^-+i)£/2, t e [-M,0], \i\ < M. 

for ah k < ko. Thus u,{t) < £{l - X* /2) with A* = ^X'^'"> and we have proved (8.61). □ 

8.4 Holder continuity 

In this section we prove the Holder regularity. 
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Theorem 8.8 Let = K with i) > 1. Suppose u is a solution to (7.57) in the time interval 
[-3A^,0] with coefficients of A satisfy (8.2), (8.3), (8.4) and (8.36) for some Z {L - K/2,L + 
K/2\ with L and K satisfying (3.1). Assume that u satisfies 

sup max|u,(t)| < £ (8.77) 
te[-ZM.a\ * 

for some £. Then there exist constants C2 > and ag > (independent of the previous exponents 
p, ^ and 1? ) such that for any < a < q;o we have 

Oscq(c)(z){u) <^X-"^", Q^"\Z) := [-Xi-",0] x {Z - M^-" , Z + M^-°'l (8.78) 

i.e. the oscillation of the solution on scale and away from the edges of the configuration 

space is smaller than tAi~'^^°'. The constants and cxq depend only on A* from Lemma 8.7 which 
depends eventually only on Eq from Lemma 8.1. 

Proof. Given A from Lemma 8.5, A* and e from Lemma 8.7, there is a small scaling factor 
V < (1/8)1/^ such that 

[{Wx\ X-% - 1]^ < i(l - \*/2)[i\x\ - X-% - 1]^ 
for any |x| > This means that 

^M/,.,z,f,A,e < 1(1 _ A72)^^f ■^'^■^■'^ if \i-Z\> M/u. (8.79) 
Set il/fc = ly'^M, and £k = (1 - A*/4)'=^, then we have 

1 ^M.,Z,£.,A,e < l^M,+i,Z,£, + i,A,e |i _ Z| > Mfe. (8.80) 



1 - A*/4 ' - 2 

Now we define the following rescaled sequence of solutions: 

u[°\t) ^u,{t), te[-3M,0] 

and 

where m^'^-' is determined by 

sup|uW - sC^'^l ^ ioscQ.u^*^) (8.81) 

(with Qk = Qm^ for brevity). 

Wc remind the reader that subtracting a constant from a solution no longer yields a solution. 
But the conclusions of previous subsections still hold as long as we made some very minor modifi- 
cations to the arguments. Since this was explained in the proof of Lemma 8.7, we will not repeat 
it here. So from now on, we will apply these Lemmas to u'"'^ . 

We claim that 

max|Mf^(0| < ik+i^f'"^^''"'^'', t e [-3Mfc,0] (8.82) 

i 

and 

OscQ.iiW < (2- A*)4 (8.83) 
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for any k. By Lemma 8.7, the first statement (8.82) implies the second (8.83). 

For fc = 0, (8.82) (hence (8.83)) certainly holds by (8.77). We proceed by induction, suppose 
(8.82) and (8.83) hold for some fc, we show (8.82) for fc + 1, i.e. we need to show that 

\^^^'\t)\ - \u^\t) - u^^^l < 4+1 t e [-3il4+i,0] (8.84) 

Using (8.83), we have 

OscQ,u('=) < (2 - A*)4 < 24+1 (8.85) 

so for \i — Z\ < Mk we have that (8.84) is satisfied by using (8.81). Moreover, we have |ti(*'')| < t^. 
Now for \i — Z\ > Mk, (8.80) can be used. Moreover, in this regime we have 



> 



-fe+1 



/ Mk 4V 



> 2^''4+i 



Choosing i/ sufficiently small {ly^ < 1/8) we can estimate this from below by Mk+i- Thus we 
estimate 



\ut\t) - ut'^)! <|«f)(t)| +4 < 24 + 



<24 + 



teh3A4+i,0], \t-Z\>Mk. 



By definition, 24 < 3£fc+i if A* is small enough. Also, we have Ut+i < if 
\i — Z\ > Mk- Thus we can bound the last equation from above by 



< 4+1 + 24+1 + 



< 4+1 + Vi 



Thus we have proved (8.82) and (8.83) by induction, i.e. 

OscQ.uW <41-A74)^- 
Applying it to a box of size Mk ^ A^^~", we have ly'^ = Ai^" and thus we obtain (8.78) with 

log(l-AV4) 



C2 



log 1/ 



This completes the proof. 



□ 



8.5 Proof of Theorem 7.7. 

With the Holder regularity result. Theorem 8.8, we now complete the proof of Theorem 7.7. Here 
we return to the notations of Section 7. 

From the decay estimate (7.42) with b = (sec (7.60)) and j|u(0)||i = 1, we have 

||u(s)||^ < s-^K^'. (8.86) 

As in the conditions of Theorem 7.7, let a G [K'^^ ^CiKXogK] be fixed and recall the relation 
M"^ K. We now choose < l/ca in Theorem 8.8 so that CM > cr > 97W (recall 7W = 2M + 1). 
Before we can apply Theorem 8.8, we have to shift the origin of the time coordinate axis by a so 
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that the time interval [— StM, 0] in Theorem 8.8 corresponds to [a — iM.,a\ in Theorem 7.7. After 
this time shift, (8.86) imphes (8.77) if we define 

e = M-^K^'. (8.87) 
From the Holder estimate, Theorem 8.8, we have 

0sCqM(2)(u) < < Ccr-l+(''«'-'=^"). gg) 

Now we choose a := minjao, \q'} and q := ^C2a, where ao was obtained in Theorem 8.8 and q' 
was given in Theorem 7.7. In particular, \j — Z\< a^~'^ implies |j — Z\ < A^^~". Then, choosing 
^0 small enough so that i?^o < q = — q, we have, since ^' < that 

Oscqm(2)(^/) <a-i-?. (8.89) 

This proves Theorem 7.7. q 

A Proof of Lemma 3.5 

First we show that on the set TZl,k, the length of the interval J = Jy = {ul-k-i, Ul+k+i) satisfies 
(3.22). We first write 

\J\ = IVL+K+l-VL-K-ll = hL+K+l-lL-K-l\+OiN-'+^'/^). (A.l) 

Then we use the Taylor expansion 

g{x) = g{y) + 0{x - y) 

around the midpoint y of J. Here we used that p G away from the edge. Thus from (2.14) 

JC + 1 = N g = N g + OiN^^/"^)^ N\J\g{y) + 0{N\J\^) + 0{N^^/^), (A.2) 

since the contribution of the second order term in the Taylor expansion is of order iV| Jp. Express- 
ing I J| from this equation and using (3.1), we arrive at (3.22). 

Now we prove (3.23). We set 

j : \j-L\>K+KS 

The potential U is similar to Vy, but the interactions with the external points near the edges of J 
(?/j's with Ij — i| < K + K^) have been removed and the external points yj away from the edges 
have been replaced by their classical value 7j. In proving (3.23), we will first compare Vy with an 
auxiliary potential U and then we compute U'. 

First we estimate the difference Vy{x) — U'{x). We fix a; G J, and for dcfinitencss, we assume 
that d{x) ^ X — t/L-K-i: x is closer to the lower endpoint of J; the other case is analogous. 
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We get (explanations will be given after the equation) 



N ^ \x — yA N ^ X — u,- X — 7," 

A'<|j-L|<A'+K« ' \j~L\>K+Ki ' 



GK^ iV-i+^«/2 1 

< 



L-K-K^ N(l-a/2) ^ ^ 



■Nd(x) d(x) N 



[ E + E 



aN/2 Ar_Ar3/f'+= 7v^/^+^ W 

[ E 1+ S E '+ E 



CK^ CjV-l+^^/^OgiV 4/15+s 

-Nd{x)^ d{x) 

Here for the first bulk sum, j e lNa/2,L - K - K% we used \yj - 7j| < iV"i+?''/2 from the 
definition of TZl,k and the fact that for j < L — K — we have 

x-jj> Vl-k-i - Ij 

> IL-K-I - Ij - \VL-K-l - lL-K-l\ 

> cN-\L -K -1- j) - CN-^+^^/^ 

> c'N-\L - K - 1 - j) 

with some positive constants c, c'. This estimate allows one to sum up \x — at the expense 

of a log factor. Similar estimate holds for j > L + K + . In the intermediate sum, j € 
lN^/^+^,Na/2\, we used \y.j - 7j| < CA^-^/is+e and that \x - y.j \ and \x - 7j | are bounded from 
below by a positive constant since 

X-yj> VL-K-l - Vj > VaN - Vj > iNa ~ lNa/2 + 0(iV"^+^*/^) > C 

and similarly for x — ^j. Finally, very near the edge, e.g. for j < N^/^+'^ ^ we just estimated \yj — 7j| 
by a constant. This explains (A. 3). 

Now we estimate U'{x). We use the fact that the equilibrium measure g = Qv satisfies the 
identity 

\v\x) = I Ady 



2 J x-y 

from the Euler-Lagrange equation of (2.13), see [1, 9]. Thus 

l\U'{x)\ < |r!i| + |r!2| + |f^3| 
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with 



L-K-K^ 



L-K-Ki g{^y) 



A 



^ L-K-K" ^ 



X — u X — 7,- 

Ti + K + ffS " j^L+x+ifS '-^ 

where i3] is the support of the density p. 
To estimate fii, we use Taylor expansion 

g{y) = gix)+0{\x-y\). 

For dcfinitencss we again assume that d{x) = x ~ Vl-k-i, and use that on TZl.k we have 

iL-K-Ki < VL-K-l <X < yL+K+l < iL+K+Ki- 

We thus obtain 

=g(x)log ^^+^^+^^'% 0(i^/jV) 

- iL^K-Ki 

=^?(y)log^ + 0(/^vv-^+«). (A.4) 

In the first step above we computed the leading term of the integral, while the other term was 
estimated trivially using that the integration length is ^i^^k+k^ ^ Il-^k-K'^ = 0{K/N). In the 
second step we used that g € away the edge, i.e. g{x) = g{y) + 0{K/N). To estimate the 
logarithm, we used 

iL+K+Ki - X ={lL+K+KS - iL+K+l) + ilL+K+l - yL+K+l) + {VL+K+l - x) 

^giy)N-^Ki + 0{N-^+^"^) + {yL+K+i - x) 
=d+{x)+0{N-^+^^/^) 

and the similar relation 

X - IL-K-K^ - d^{x) + 0(iV-i+«^/2). 

Notice that the error term in (A.4) is smaller than the target estimate / {N d{x)) since d{x) < 
K/N < K-^+^N-^. 

Now we estimate the ^2 term; fis can be treated analogously. We can write (with the convention 
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70 = 



73 



E 



£'(y) 



1 



1 



<C E (7,-7,-i) r T 

j = l •'^3-1 



dy 



<CN-'^ 



dy 



A + K 



{x - yY 



X - y) 



A+K 



dy 



{x - yy 



<- 



CN- 



d{x) 

In the first step we used that 



7j 



73-1 



1 

N 



from (2.14). In the second step we used that 7j — 7j_i = Ok{N~^) in the bulk, i.e. for > A + k, 
and niaxj(7j — 7j_i) = 0(iV~^/^) (the order A^"^/^ comes from the fact that the density p vanishes 
as a square root at the endpoints). The parameter k = K{a) is chosen such that A + 2k < yL-K-i 
which can be achieved since L > aN and yL-K-i is close to 7l-x-i- In the very last step we 
absorbed the error term into {Nd{x))~^ > K^^ > N~'^/^. 

Finally we prove (3.24). Since \yj — jj\ < K^/N, it follows that \x — yj\ ~ | 
|a; — 7il > K^/N. Thus we have 



for 



" Nf^^ix-y^y- 7V^(x-7j)2 



> inf V" 



d{x) 



with some positive constant c (depending only on a). In estimating the summation, we used that 
the sequence 7^ is regularly spaced with gaps of order 1/7V. This completes the proof of Lemma 3.5. 

□ 

B Critical discrete Gagliardo-Nirenberg inequality for \p\ 

Proposition B.l There exists a positive constant C such that 



\n\Uz)<c\\f\\u 



J2 l/^:-/.f 



;|2 



(B.l) 



holds for any function / : Z — > M. 

Proof. Recall the integral formula for quadratic form of the operator \p\ in R [39] 

1 f f \cf>{x)-c^{yr 



(^(x) {\p\ (j)){x)dx 



27r 



\x - 



-dxdy. 



(B.2) 
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We have the the following Gagliardo-Nirenberg type inequality in the critical case (see (1.4) of [42] 
with the choice n = l,p = A) 

\mt<C\ml f ct>{x){\p\m^)dx, cP-.R^R. (B.3) 
Jr 

Thus (B.l) is a discrete version of (B.3). Now we show how to to derive (B.l) from (B.3). 
Given / : Z M, let : M M be its linear interpolation, i.e. := fi for i G Z and 

Hx) = /. + (/.+! - -t) = h+i - - h){i + 1 - x), x€[i,i + l]. (B.4) 

It is easy to see that 

C^~'I10||lp(k) < <C||0|Up(r) , 2<p<4, (B.5) 

with some universal constant C . We claim that 

with a universal constant C, then (B.5) and (B.6) will yield (B.l) from (B.3). 
To prove (B.6), wc can write 

Using the explicit formula (B.4), we first compute the i = j terms in (B.7): 

Next we compute the terms |i — j| ~ 1 in (B.7). We assume j =2 — 1, the terms j = i + 1 are 
analogous; 

where we used (l){x) = fi + (Z^+i - fi){x - i) and 0(j/) = - (/; - fi-i){i - y)- The above integrals 
are finite constants, so we get 



i-l 



\x-y\'^ " - {i + l-if {i-{i-l)f 



Finally, for the terms |i — j| > 2, we can just replace (a; — y)'^ by {i — jf' in the right hand side of 
(B.7) and use simple Schwarz inequalities to get 



\x-y\^ 



2 
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After summing up we get 

The estimates (B.8), (B.9) and (B.IO) together yield (B.6). □ 
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